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ABSTRACT 


Physical parameters of both the mass donor and compact object can be constrained in X-ray 
binaries with well-defined eclipses, as our survey of wind-fed supergiant X-ray binaries (SGXBs) 
IGR J16393-4643, IGR J16418-4532, IGR J16479-4514, IGR J18027-2016 and XTE J1855-026 
reveals. Using the orbital period and Kepler’s third law, we express the eclipse half-angle in 
terms of radius, inclination angle and the sum of the masses. Pulse-timing and radial velocity 
curves can give masses of both the donor and compact object as in the case of the “double-lined” 
binaries IGR J18027-2016 and XTE J1855-026. The eclipse half angles are 15^^, 31.7t5^-^, 32±2, 
34±2 and 33.6±0.7 degrees for IGR J16393-4643, IGR JI64I8-4532, IGR JI6479-45I4, IGR 
J18027-2016 and XTE 1855-026, respectively. In wind-fed systems, the primary not exceeding 
the Roche-lobe size provides an upper limit on system parameters. In IGR J16393-4643, spectral 
types of BO V or BO-5 III are found to be consistent with the eclipse duration and Roche-lobe, 
but the previously proposed donor stars in IGR J16418-4532 and IGR J16479-4514 were found 
to be inconsistent with the Roche-lobe size. Stars with spectral types 07.5 I and earlier are 
possible. Eor IGR J18027-2016, the mass and radius of the donor star lie between 18.6-19.4 Mq 
and 17.4-19.5i? 0 . We constrain the neutron star mass between 1.37-1.43 Mq. We find the mass 
and radius of the donor star in XTE J1855-026 to lie between 19.6-2O.2M0 and 21.5-23.Oi?©. 
The neutron star mass was constrained to 1.77-1.82M 0 . Eclipse profiles are asymmetric in IGR 
J18027-2016 and XTE J1855-026, which we attribute to accretion wakes. 


1. Introduction 

High-Mass X-ray Binaries (HMXBs) are rel¬ 
atively young systems, which consist of a com¬ 
pact object (neutron star or black hole) and an 
early-type OB star orbiting the common center 
of mass. Eirst discovered in the 1970s, HMXBs 
are split into two individual classes-Be X-ray bi¬ 
naries (BeXBs) and supergiant X-ray binaries 
(SGXBs). BeXBs are transient systems where 
the compact object is in a wide, typically eccen- 
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trie, orbit (P>10days) around a rapidly rotating 
non-supergiant B-type star. The primary mode 
of mass transfer occurs when the compact object 
passes through the circumstellar decretion disc 
of the mass donor. In SGXBs, the compact ob¬ 
ject is in a short (~l-42day) orbit around an 
OB-supergiant where the accretion mechanism 
is either via the powerful stellar wind and/or 
Roche-lobe overflow. While the eccentricity in 
most SGXBs with short orbital periods is near 
zero, some SGXBs host compact ob jects in highly 
eccen tric orbits (e.g. GX 301-2 llslam fc Pauli 
20141) . In wind-accretors, the X-ray luminosity 
is on the order of 10^®-10^®erg s“^. However, 
a much higher luminosity is found in systems 
where the donor fills its Roche-lobe, ~10^®erg 


s ^ (iKaoer et al. 


Many HMXBs host a neutron star where mod¬ 
ulation is often seen at its rotation period. The 
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mass of the neutron star can be constrained in 
eclipsing X-ray pulsars, which can lead to an im¬ 
proved unders tanding of the neu tron star Equa¬ 


tion of State (|Mason et alJ l201lL and references 


therein). Currently, the neutron star mass has 
been constrained in 10 XRBs where the lower limit 
for each object corresponds to edge-on orbital in¬ 
clinations and the upper limit is calculated when 
the donor star just fills its Roche-lobe. While 
over 100 Equatioii s of State have been proposed 
( Kaoer et al.l 2006ll. only one mo del is physically 
correct ( Mason et al. 20 id 20111 . 

The mass ratio between the neutron star and 
donor star is equal to the ratio between the semi¬ 
amplitude of the radial velocities of the donor 
star, Kq, and the neutron star, Kx- Through¬ 
out the paper, we use a definition of mass ra¬ 
tio, q, as the ratio of the co mpact object mass 
to th at of the donor star mass (j.Toss fc Rapnaport 


198411 . The orbital period of the binary, P, and 


Kx can be calculated us ing pulse-timing analy¬ 
sis ( Val Baker et al.ir2005 , and references therein), 
which is analogous to measuring the Doppler shift 
of spectral lines in the opt ical and/or near-infrared 
( Joss fc Rappaporll Il984 1. The projected semi¬ 
major axis can be determined from the semi¬ 
amplitude of the radial velocity of the neutron 
star. The semi-amplitude of the radial velocity 
of the donor star can be determined using optical 
and/or near-infrared spectroscopic information. 

Eclipse measurements can also be exploited as 
timing markers to determine the binary orbital 
evolution of HMXBs. A significant orbital pe¬ 
riod derivative, P, was previously found in several 
eclipsing HMXBs (e.g 4U 1700-377; S MCX-1; Gen 


X-3; LMC X-4 and OAO 1657-415. iRubin et al 


Il996l : iRaichur fc Paul 2010t iFalanga et al.l l2015ll 
and can be used to investigate the orbital evo¬ 
lution over long periods of time. While several 
contending theories to explain the orbital decay 
have been investigated, the most probable expla¬ 
nations involve tidal interaction and rapid mass 
transf er between the comp onents of the binary sys¬ 
tems (jEalanga et al.ll2015L and references therein). 

The Swift Burst Alert Telescope ( BAT), sensi¬ 


tive t o X-rays in the 15-150 keV band (|Barthelmv et al. 


20051) . provides an excellent way to study highly 


absorbed SGXBs. The large absorption found in 
these systems is problematic for instruments such 
as the Rossi X-ray Timing Explorer (RXTE) All 


Sky Monitor (ASM), which ope rated in the 1.5- 
12keV band (ILevine et al.lPOllI ). The sensitivity 
to higher energy X-rays a llows Swift-BAT to peer 
through this absorption (ICorbet &: KrimmI 12013 , 
and references therein). 

We present here constraints on the mass and 
radius of the donor star in eclipsing XRBs. The 
probability of an eclipse in a supergiant XRB with 
an orbital period less than 20 d can be expressed 
in terms of the orbital period, mass of the donor 
star, and radius of the d onor star (Equation 1; 
Bay fc Chakrabartv 20021) . The probability of an 
eclipse in long-period XRBs is low. Using a litera¬ 
ture search, we determined HXMB systems where 
BAT observations can significantly improve the 
the properties of the stellar components and or¬ 
bital evolution of the systems using eclipsing prop¬ 
erties. Five eclipsing XRBs were identified: IGR 
J16393-4643, IGR J16418-4532, IGR J16479-4514, 
IGR J18027-2016 and XTE J1855-026, which are 
all highly obscured SGXBs. We note while the 
masses of both the donor and compact object had 
previously been constrained in IGR J18027-2016, 
the error on the eclipse half-angle was large at 4.5° 
( Hill et al]l2005l) . The error estimates concerning 
the radius and mass of the donor star are signifi¬ 
cantly improved in our analysis. 

This paper is structured in the following order: 
Swift BAT observations and the description of the 
eclipse model are presented in Section 2; Section 
3 focuses on individual sources that are known to 
be eclipsing. Section 4 presents a discussion of the 
results and the conclusions are outlined in Section 
5. If not stated otherwise, the uncertainties and 
limits presented in the paper are at the Itr confi¬ 
dence level. 

2. Data Analysis and Modelling 
2.1. Swift BAT 

The BAT on board the Swift spacecraft is 
a hard X-ray telescop e operating in the 15 - 
150keV energy band ( Barthelmv et al.1 20051) . 
The detector is composed of CdZnTe where 
the detecting area and field of view (FOV) are 
5240cm^ and 1 . 4sr (half-coded), respectively 
( Barthelmv et al. 20051) . The BAT provides an 


all-sky har d X-ray survey wit h a sensitivity of 


^ImCrab (jTueller et al 
duces ~0.045 counts cm" 


201fll) . The Grab pro- 


^ s ^ over the entire en- 
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ergy band. 

We analyzed BAT data obtained during the 
time period MJD 53416-56745 (2005 February 15- 
2014 March 29). Light curves were retrieved using 
the extraction of the BAT transient monitor data 
available on the NAS A GSFC HEASARC website0 
( Krimm et al. |[20 i 3) , which includes orbital and 
daily-averaged light curves. We used the orbital 
light curves in the 15-50 keV energy band in our 
analysis, which have typical exposures of ^6 min 
(see Section [3]). The short exposures, which are 
somewhat less than typical Swift pointing times 
(^20min), can arise due to the observing plan of 
Swift itself where BAT i s primarily tasked t o ob¬ 
serve gamma-ray bursts ( Krimm et al.ll^OlSl) . 


The light curves were further screened to ex¬ 
clude bad quality points. We only considered data 
where the data quality flag (“DATA_FLAG”) was 
set to O. Data flagged as “good” are sometimes 
suspect, where a small number of data points with 


verv low fluxes and implausibb 

i small uncertainties 

were found ( 

Corbet & Krimm 

2013). These Doints 


were removed from the light curves. We cor¬ 
rected the photon arrival times to the solar system 
barycenter. We used the scripts made available on 
the Ohio State Astronomy webpageH In this pa¬ 
per, the barycenter-corrected times are referred to 
as Barycenter Modified Julian Date (BMJD). 


We initially derived the orbital period for each 
XRB in our sample using Discrete Fourier trans¬ 
forms (DFTs) of the light curves to search for 
periodicities in the data. We weighted the contri¬ 
bution of each data point to the power spectrum 
by its unc ertainty, using the “semi-weighting” 


technique (jGorbet et al.l l2007t ICorbet fc Krimm 


2013ll . where the error bars on each data point 


and the excess variability of the light curve are 
taken into account. We derived uncertainties on 
the orbital period s using the expression given in 
Horne fc Balinna^ (1986). 


2.2. Eclipse Modeling 

We initially modeled the eclipses using a sym¬ 
metric “step and ramp” function (see Table [T]) 
where the intensities are assumed to remain con¬ 
stant before ingress, during eclipse and after egress 
and follow a linear trend during the ingress and 

^ http: //heasarc.gsfc.nasa.gov/docs / swift / results / transients / 
^ http: / / astroutils. astronomy, ohio-state.edu / time /. 


egress transitions ( Colev et al.ll2014ll . The count 
rates before ingress, during eclipse and after egress 
were each considered to be free parameters and 
were fit as follows: Cing was ht from binary phase 
(j) =-0.2 to the start of ingress, Ced was fit during 
eclipse and Ceg was fit from the end of egress to 
phase (j) =0.2 (see Equation[T]). While we find the 
eclipse profiles of IGR J16393-4643, IGR J16418- 
4532 and IGR J16479-4514to be symmetric within 
errors, the profiles show some asymmetry in the 
cases of IGR J18027-2016 and XTE J1855-026. 


We note that a symmetric “step and ramp” func¬ 
tion could lead to systematic errors and we there¬ 
fore fit the eclipse profiles using an asymmetric 
“step and ramp” function (see Equantion [T] and 
Table [2]). The parameters in this model are as 
follows: the phases corresponding to the start of 
ingress and start of egress, ^ing and ^eg, the du¬ 
ration of ingress, A(/)ing, the duration of egress, 
A())eg, the pre-ingress count rate, Cing, the post¬ 
egress count rate, Ceg, and the count rate during 
eclipse, Ced. The eclipse duration and mid-eclipse 
phase are calculated using Equations [5] and [31 re¬ 
spectively. 


Cing, 



-0.1 

^ < (/> < 4>ing 


^ing) ^ingj 

(l>ing 

< 0 < 0ing + A(/iing 

Cecl, 



(l>ing 

+ A0ing ^ (p ^ 0egr 


0egr) Cecil 

4^egr 

^ ^ ^ 0egr A0eg 

Ceg, 



0egr 

Al^eg ^ (j) 0.2 





(1) 


— ^egr 

— (<('ing + A(/iing) 

(2) 

4>mid = 

2 ('/’egr 

+ (</'ing + A(/)ing)) 

(3) 


The eclipse duration, time of mid-eclipse, and 
eclipse half-angle (0e = A())edXl80°) from htting 
the BAT folded light curves for each source are 
reported in Tables [T]-[31 For each source, we ini¬ 
tially used an ephemeris based on our determi¬ 
nation of the orbital period from the DFT and 
time of mid-eclipse. Using an ‘observed minus 
calculated’ O — C analysis (see Figures |TH3]), we 
refined the orbital periods and improved on the 
time of mid-eclipse for each XRB in our sample. 
We note that no eclipses are visible in the un¬ 
folded light curves and it is necessary to observe 
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Table 1 

Eclipse Model Parameters, Assuming a Symmetric Eclipse Profile 


Model Parameter 

IGR J16393-4643 

IGR J16418-4532 

IGR J16479-4514 

IGR J18027-2016 

XTE J1855-026 

4’ing 

-0.0791H?! 

-0.107T0.002 

-0 1 14+0-003 
'^•-*--*-^-0.004 

U.10U_Q 003 

-0 141+0-001 

U.J.OJ._o 002 

A4> 

p 

o 

o 

1 + 
o o 
b b 
o o 

00 CD 

0 01 Q+0 002 

u.uiy_o.oo3 

0 noQ+o oos 

u.uzy_o.oo4 

U.Uc)0_o.oo3 

0.038T0.002 

(7“ 

1.27T0.04 

1.18±0.05 

1.05±0.06 

1.59T0.06 

2.64±0.06 

(pegi 

^ iD 
O O 

o q 
b b 
+ 1 
00 

o 

d 

0.088l°°°? 

0 nqo+o oos 
u.uyz_Q QQ4 

O.OOSlHoI 

0.094±0.002 

a 

^ecl 

0.71T0.06 

0.00T0.05 

0.03±0.05 

0.17T0.04 

-0.04T0.05 

^0ecl 

o.o97;;:;^2 

0.175±0.003 

n 1 

U.i M _0.007 

(TTo^rnTruT 

U.iyD_o.o05 

0.187T0.003 

pb 

orb 

4.23794±0.00007 

3.73880T0.00002 

3.31965T0.00006 

4.56999T0.00005 

6.07410T0.00004 

pc 

-5±4 

0.7T1.0 

3±2 

-2±2 

0.5±1.0 

rj-id 

mid 

55074.99l|^'24 

55087.714T0.006 

55081.571 

55083.79tg);i 

55079.055t);j(o^ 



31.5T0.6 

3i-9l?:3 

35±1 

33.61°:® 

Xt (dof) 

1.13(77) 

0.84(77) 

1.03(77) 

0.93(77) 

1.23(77) 


Note.— 

“ Units are 10“® counts cm~^ s“^. 

^ Refined orbital periods using an O — C analysis. Units are days. 

The orbital period derivative at the 90% confidence interval found using an O — C analysis. Units are 10~^d d~^ 
Units are barycentered Modified Julian Time (BMJD). Phase 0 is dehned as eclipse center. 

® Units are degrees. 


Table 2 

Eclipse Model Parameters, Assuming an Asymmetric Eclipse Profile 


Model Parameter 

IGR J16393-4643 

IGR J16418-4532 

IGR J16479-4514 

IGR J18027-2016 

XTE J1855-026 

4’ing 

-0.092lg;°o5 

-0.108T0.002 

-0.118T0.003 

-0.160T0.005 

-0.136T0.002 

^0ing 

0.06T0.01 

0.020T0.003 

0.024T0.008 

0.08±0.01 

0.042T0.003 

A((>eg 

0.039T0.006 

0 014+0-003 
U.Uilt_0.004 

U.UZD_0.003 

0.018ll;“i 

0.038T0.003 

^ a 

Oing 

1.28±0.06 

1 07+0-07 

J-.U^_0.06 

1.01T0.08 

1.43T0.08 

2.44±0.08 

a 

Oeg 

1.26T0.06 

1.22±0.07 

1.04±0.08 

1.67T0.08 

Q '7Q + 0.09 

1 ^-0.08 

4>egi 

u.uity_o.oo6 

0.089T0.002 

0.086T0.003 

0.112±0.002 

u.uyz_Q QQ2 

a 

^ecl 

0.69T0.06 

0.00T0.05 

0.05±0.05 

0.14T0.05 

-0.05±0.05 

^0ecl 

0.09Tl,^t 

0.176T0.004 

0.180±0.009 

0.19T0.01 

0.186±0.004 

pb 
^ orb 

4.23810T0.00007 

3.73881T0.00002 

3.31961T0.00004 

4.56988±0.00006 

6.07413±0.00004 

pc 
^ orb 

-6±4 

0.2T0.8 

3±2 

-3±4 

0.5±1.9 

rj-id 

mid 

55074.99T0.03 

55087.72111;“^ 

55081.57T0.02 

55083.88±0.03 

55079.07T0.01 

es 

1511 

31 - 711 ;^ 

32±2 

3411 

33.6T0.7 

xi (dof) 

1.18(80) 

0.97(80) 

1.08(80) 

1.01(80) 

1.03(80) 


Note.— 

“ Units are 10“® counts cm~^ s“^. 

^ Refined orbital periods using an O — C analysis. Units are days. 

The orbital period derivative at the 90% confidence interval found using an O — C analysis. Units are 10~^d d~^ 
Units are BMJD. Phase 0 is dehned as eclipse center. 

® Units are degrees. 
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Table 3 

Eclipse Model Parameters, Including Historical Mid-Eclipse Times 


Model Parameter 

IGR J18027-2016“'^ 

IGR 318027-2016“’^ 

XTE 31855-026*'’“ 

XTE 31855-026*'’'* 

0ing 

-0 147+0004 
'4-4+' -0.005 

0 1 A7+0.005 
U.iO/_0.004 

-0 1*41+0001 
U.lOl_0.002 

-0.136±0.002 

A(/)ing 

U.UOO_Q QQ7 

0.082±0.009 

0.038±0.002 

0.040±0.003 

A'/'eg 

U.UOO_Q QQ7 

U.UZ / _0.006 

0.038±0.002 

0.037±0.003 

r'. e 

'-^ing 

1.59±0.05 

1.25±0.08 

2.64±0.06 

2.45±0.08 

Oeg 

1.59±0.05 

1.68±0.08 

2.64±0.06 

9 7Q+0.09 
'1^-0.08 

4^egr 

0.099±0.003 

0.102±0.002 

0.094±0.002 

u.uyz_o 002 

Cecf 

0.17±0.04 

0.14±0.04 

-0.04±0.05 

-0.03±0.05 

^^ecl 

U.iyO_Q QQg 

0.19±0.01 

0.187±0.003 

0.187±0.004 

pf 

4.56982±0.00003 

4.56993±0.00003 

6.07412±0.00003 

6.07414±0.00003 

p9 

-^orb 

0.8±0.9 

0.2±1.1 

-0.1±0.5 

0.0±0.5 

rph 

55083.78±0.01 

55083.82±0.01 

55079.056±0.009 

55079.07±0.01 


35±1 

34±2 

33.6t°j 

33.7±0.7 

xH^) 

1.20(77) 

0.93(80) 

1.02(77) 

1.05(80) 


Note.— 


(12005 ). Jain et al. ( 2009) and Falanga et al. ( 2015 ). 


“ Includes the mid-eclipse times derived in iHill et al.l 

^ Includes the mid-eclipse times derived in Icorbet fc Mukail ( 2002 ) and Falanga et al. ( 20151) . 

Assuming a Symmetric Eclipse Profile. 

Assuming an Asymmetric Eclipse Profile. 

® Units are 10“^ counts cm“^ s“^. 

f Refined orbital periods using an O — C analysis. Units are days. 

® The orbital period derivative at the 90% confidence interval found using an O — C analysis. Units are 
10“^d d“^. 

^ Units are BMJD. Phase 0 is defined as eclipse center. 

* Units are degrees. 
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multiple cycles of folded light curves in order for 
eclipses to be seen. We divide the light curves 
into five equal time intervals (~670 days), with the 
exception of IGR J16393-4643 (see Section IXTl) . 
and calculate the mid-eclipse epoch for each in¬ 
terval (see Table lU In the cases of IGR JI8027- 
2016 and XTE J1855-026, we combine our derived 
mid- eclipse times with tho se reported in the litera¬ 
ture (IFalanga et al.ll2015l . and references therein). 
We note that while mid-eclipse times were previ¬ 
ously derived for both IGR J16393-4643 and IGR 
J16479-4514 (see Table E]), these were not used 


since no error estimate was reported (jlslam et al 


2 OI 5 I : Bozzo et al. 2009ll . We then fit the mid¬ 
eclipse times using the orbital change function (see 
Equation |4]) where n is the number of binary or¬ 
bits given to the nearest integer, Porb is the orbital 
period in days, Porb is the period derivative at Pqj 
and the error on the linear term is the orbital pe¬ 
riod error. In all five cases, we improve the error 
estimate on the orbital period by nearly an order 
of magnitude (see Section [3]). We do not find a 
significant Porb for any source (see Tables [U-El). 


Tn — Tq nPorb + TT^^^-forbT^orb 


( 4 ) 


X-ray binaries that are eclipsing have an eclipse 
duration that is only dependent on the radius of 
the mass donor, inclination angle of the system 
and the orbital separation of the components pro¬ 
vided that the orbit is circular (e = 0). Us¬ 
ing the observed orbital period and Kepler’s third 
law, the duration can be written in terms of the 
sum of the donor star and compact object masses, 
which stipulates that the eclipse half-angle, 0e, 
can now be expressed in terms of the radius, in¬ 
clination and masses of the components. In one 
set of calculations, we assume a 1.4 M© compact 
object which ma y be appropriate for an accret¬ 
ing neutron star (jGhandrasekharlll931r i. The re¬ 
gion allowed by the measured eclipse half-angle 
for each binary in Mass-Radius space is shown 
in Section [ 3 ] The inclination is constrained be¬ 
tween edge-on orbits (left boundary of the dark 
shaded region) and close to face-on orbits (the 
right boundary of the light shaded region). We 
can attach additional constraints assuming that 
the mass donor underfills the Roche-lobe (right 
boundary of the dark shaded region), which is de¬ 
pendent on the mass ratio of the system and the 


u 

I 

o 

"cd 

T3 


0.0 

-0.1 

0.02 

0.00 

-0.02 

0.1 

0.0 

0.10 

0.05 

0.00 

-0.05 

0.04 

0.02 

0.00 

- 0.02 


IGRJ16393-4643 


1 


IGRJ16418-4532 


1GRJ16479-4514 


1GRJ18027-2016 - 


XTEJ1855-026 - 


H- 


-400 -200 0 200 400 

Cycle 


Fig. 1.— The observed minus calculated (O — C) 
eclipse time residuals for IGR JI6393-4643 (top), 
IGR J16418-4532 (second panel), IGR JI6479- 
4514 (middle), IGR J18027-2016 (fourth panel) 
and XTE J1855-026 (bottom) fit using a symmet¬ 
ric step-and-ramp function. We subtract the best 
linear polynomial fit for each source and correct 
the orbital periods accordingly. For IGR J16393- 
4643 (top) we only use the first four points to ob¬ 
tain a good fit (see Section [3d1) . 


orbital separation. To calculate the eclipse half¬ 
angle an d the Roche-lobe radius, we used Equa- 
tion 7 in Joss fc Rappaport ( 1984ll . als o used by 


Eggleton 


iRubin et al.l (Il996ll . and ^uation 2 
(|1983ll . respectiyely. Further constraints on the 
parameters of the donor star are imposed with 
pulse-timing techniques (dashed red lines in Fig¬ 
ures and [HI). For the systems where pulse¬ 
timing results were not ayailable, we additionally 
calculated the minimum inclination angle of the 
system, imin, that is consistent with the measured 
eclipse half-angle (see Tabled]). 

When the semi-amplitude of the radial yeloc- 
ities of both the compact object and the mass 
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Table 4 

Mid-eclipse Time Measurements For O - C Analysis 


Source 

Orbital 

Mid-eclipse time“ 

Mid-eclipse time*” 


Cycle (N) 

(MJD) 

(MJD) 

IGR J16393-4643 

-313 

53748.46T0.03 

53748.44T0.03 

IGR J16393-4643 

-155 

54418.15in2 

54418.141°;°^ 

IGR J16393-4643 

1 

55079.22+°g^ 

55079.22lg-® 

IGR J16393-4643 

158 

55744.55l°;°2 

55744.54l°;°i 

IGR J16393-4643 

315 

C 

C 

IGR J16418-4532 

-359 

53745.50T0.01 

53745.49lll^ig 

IGR J16418-4532 

-179 

54418.4611^-^?^ 

54418.468±0.009 

IGR J16418-4532 

0 

55087.71±0.01 

55087.72T0.01 

IGR J16418-4532 

179 

55756.963l?-?n« 

55756.967±0.009 

IGR J16418-4532 

357 

56422.47±-0.02 

56422.47lK;2? 

IGR J16479-4514 

-402 

53747.191^ 

53747.211“'^ 

IGR J16479-4514 

-200 

54417.63±0.02 

54417.611°-“ 

IGR J16479-4514 

0 

55081.47l2-!(^ 

55081.48i;(-23 

IGR J16479-4514 

200 

55745.331^4 

55745.28±0.04 

IGR J16479-4514 

401 

56412.60lo;o? 

56412.59lHi 

IGR J18027-2016 

-292 

53749.411^ 

53749.401^-^5 

IGR J18027-2016 

-143 

54430.32T0.02 

54430.34l!(-g^ 

IGR J18027-2016 

0 

55083.791^3 

55083.851003 

IGR J18027-2016 

146 

55751.02T0.03 

55751.06looi 

IGR J18027-2016 

292 

56418.12l°-°2 

56418.201°;“ 

XTE J1855-026 

-219 

53748.82T0.01 

53748.84±0.01 

XTE J1855-026 

-109 

54416.99T0.02 

54416.99T0.02 

XTE J1855-026 

0 

55079.05±0.01 

55079.09l!(-23 

XTE J1855-026 

110 

55747.20T0.01 

55747.21T0.02 

XTE J1855-026 

219 

56409.30T0.02 

56409.33T0.02 


Note.— “ Obtained using the symmetric step-and-ramp function. 

** Obtained using the asymmetric step-and-ramp function. 

A bad fit for IGR J16393-4643 was obtained in the mid-eclipse time between 
MJD 56079-56745. 
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Table 5 

Historic Mid-eclipse Time Measurements For IGR J18027-2016 and XTE J1855-026 


Source 

Orbital 

Mid-eclipse time 

Satellite 

Reference 


Gycle (N) 

(MJD) 

(MJD) 


IGR J16393-4643 

-391 

53417.955 

Swift 

Islam et al. (20151 

IGR J16479-4514 

-161 

54547.05418 

Swift 

Bozzo et al. (2009) 


IGR J18027-2016 

-638 

52168.22±0.12 

BeppoSAX 

IGR J18027-2016 

-638 

52168.26T0.04 

BeppoSAX 

IGR J18027-2016 

-471 

52931.37±0.04 

INTEGRAL 

IGR J18027-2016 

-399 

53260.37T0.07 

INTEGRAL 

IGR J18027-2016 

-286 

53776.82±0.07 

Swift 

IGR J18027-2016 

-267 

53863.10±0.14 

INTEGRAL 

IGR J18027-2016 

-127 

54503.38±0.07 

SwiR 

XTE J1855-026 

-590 

51495.25±0.02 

RXTE 

XTE J1855-026 

-391 

52704.04±0.05 

INTEGRAL 

XTE J1855-026 

-176 

54009.97±0.05 

INTEGRAL 

XTE J1855-026 

-31 

54890.68±0.05 

INTEGRAL 


MigelJoe^^ (2003) 
(2005) 
(2005) 
^ain_e£_^ (2009) 
^ain_e£_^ (2009) 
(2015) 

^ain e£^ (2009) 


^orbet_&_Mukai (2002) 
JFalAiiga^_e£_^ (2015) 
^]alAn^a^_e£_^ (2015) 
FalAn^a^_e£_^ (2015) 


Note.— Historical mid-eclipse times for IGR J16479-4514, IGR J18027-2016 and XTE J1855- 
026 found using RXTE, Swift BAT and INTEGRAL. 


donor are known (e.g. IGR JI8027-20I6, XTE 
J1855-026), the mass ratio between the compact 
object a nd mass donor can be ca lculated (Equa¬ 
tion 6, IJoss fc Rappaport 1984 ). This means 
that in addition to the radius and mass of the 
donor star, the mass of the compact object can 
be constrained. The mass of the donor star 
can be written in terms of the semi-amplitude 
of the radial velocity of the compact object, or¬ 
bital period, Newton’s gravitational constant, in- 
clination angle of the sy stem and the mass ratio 
( Joss fc RapDaport 1984h . Likewise, the compact 


object mass can be written in terms of the semi¬ 
amplitude of the of the radial velocity of the donor 
star, orbital period, Newton’s gravitational con¬ 
stant, inclination angle of the sy stem and the mass 
ratio ( Joss fc Rappaport 1984f) . To calculate the 
masses of both the donor star and th e compact 


object , we used Eq uations 2 and 3 in I Ash et al 


( 1999 I) . also used bv iRappaport &: Josd ( 19831) . 


clination angle using the mass and radius for the 
derived spectral type s (see Section [3]). Gene rally 
we used results from Carroll fc Ostliel ( 2006ll for 
main-sequence, giant and supergiant luminosity 
classes. These are represented by the red, green 
and blue dashed lines in mass-radius space for each 
system (see Section[3]) . Fo r 0-type supergiants, we 
also use Tables 3 and 6 in iMartins et al.l (l2005l) to 


compare our results (see blue dotted lines in Fig¬ 
ures EM and ini). The constraints for B-type su¬ 
pergiants are additionally com pared with Tables 
3 and 6 in iLefeyer et al.l ( 2007l l for B-type super¬ 
giants (blue crosses in Figures [15] and [18]). 

3. Five Eclipsing HMXBs 

3.1. IGR J16393-4643 (=AX J16390.4- 
4642) 

IGR J16393-4643 is an HMXB first discovered 
and listed as AX J16390.4-4642 in the ASCA Faint 


For consistency, we compare our derived con- 

Source Catalog ( 

Sugizaki et al. 

2001) and was 

straints on the masses and radii of the donor 

later detected with INTEGRAL 

Bodaghee et al. 


stars with those expected for the previously pro¬ 
posed spectral types. For the systems where pulse¬ 
timing results were not available, we calculate the 
predicted eclipse half-angles as a function of in- 


The average flux in the 20-40 keV band 

and 


200611 . 

was found to be 5.1x10“^^ erg cm“^ s ", ana in¬ 
tensity variations were found to exceed a factor 
of 20 ( Bodaghee et al. 2006). In the 2-10 keV 
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Fig. 2.— The observed minus calculated (O — C) 
eclipse time residuals for IGR J16393-4643 (top), 
IGR J16418-4532 (second panel), IGR J16479- 
4514 (middle), IGR J18027-2016 (fourth panel) 
and XTE J1855-026 (bottom) ht using an asym¬ 
metric step-and-ramp function. We subtract the 
best linear polynomial fit for each source and cor¬ 
rect the orbital periods accordingly. For IGR 
J16393-4643 (top) we only use the first four points 
to obtain a good fit (see Section ITT]) . 


unabsorbed flux was found 


energy band, the 

to be 9.2x10“^^ erg cm“^ s“^ (iBodag-hee et al 


2006l l. A proposed mass donor 2MASS J16390535- 
4242137 was found in the XMM-Newton error 
circle, which is though t to be an OB supergiant 


(jBodaghee et al.M2006ll . However, a precise po¬ 
sition of the donor star obtained with Chandra 
shows this candidate to b e positionally inconsis - 
tent with the X-ray source ( Bodaghee et al.ll2012 b 
Using the Spitzer Galacti c Legacy Infrared Mid- 
Plane Survey (GLIMPSEl. lBodaghee et ^ ( 2012 ) 
proposed that the counterpart must be a distant 
reddened B-type main-sequence star. 


Fig. 3.— Derived O — C eclipse time residuals for 
IGR J16393-4643 (top), IGR J16418-4532 (second 
panel), IGR J16479-4514 (middle), IGR J18027- 
2016 (fourth panel) and XTE J1855-026 (bottom) 
fit combined with those in the literature. We sub¬ 
tract the best linear polynomial fit for each source 
and correct the orbital periods accordingly. For 
IGR J16393-4643 (top) we only use the first four 
points to obtain a good fit (see Section [XT]) . 


( 20061) found a 912.0±0.1s modulation, which 
was int erpreted as t he ne utron star rotation pe¬ 


riod. Ilslam et al.l ([2015|) recently refined this 
to 908.79±0.01s using Suzaku. A ^3.7 day or¬ 
bital period was suggested using a pulse tim¬ 
ing analysis, although orbital periods of ^50.2 
and ^8.1 days were n ot completely ruled out 


( Thompson et al.l l2006l) . While various possible 
orbital solutions and accretion mechanisms have 
been proposed, orbital periods of 4.2368±0.0007 
and 4.2371±0.0007days were clearly found from 
data from Swift-BAT and R XTE PGA, re¬ 
spectively ( Gorbet et al.l 201(il) 


This was re¬ 


using INTEGRAL and XMM-JVewton. iBodaghee et al 


hned to 4.2386±0.0003d ( Gorbet fc Krimn]l2013 ). 
also using BAT data. Ilslam et ahl ( 2015h re- 
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cently derived an orbital period of ~366150s 
(4.24 d) with B AT, which is also c o nsiste nt with 
the result from Corbet fc Krimnil ( 201,3 1. The 
position in Corbet’s di agram shows that IGR 
J16393-4643is anSGXB ( Corbet fc Krimmll2013 1. 
Corbet fc Krimnil ( 2013 ) identified the presence of 
a possible superorbital period of '^ISdays; al¬ 
though with low significance. 
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Fig. 4.— Swift-BAT light curve of IGR J16393- 
4643 in the 15-50keV band folded on the or¬ 
bital period (top) using 20 bins. TO is defined 
at BMJD 55074.99, corresponding to mid-eclipse. 
A detailed folded light curve with 80 bins (bot¬ 
tom) is fit with both a symmetric “step and ramp” 
function (green) and asymmetric “step and ramp” 
function (red), which model the eclipse. The sym¬ 
metric “step and ramp” function was shifted ac¬ 
cordingly. 


The BAT light curves folded on the orbital pe¬ 
riod revealed a sharp dip, which was interpreted 


as an eclipse (jCorbet fc Krimmll2013t llslam et al 
2 OI 5 II . With Swift BAT, Islam et al. ( 2015ll cm- 
strained the eclipse half angle to be ^17°, cor¬ 
responding to a duration of ^0.75d (~65.1ks). 
Using the relationship between eclipse duration 
and stellar radius, a long with the defin i tion o f 
the Roche-l o be fr om Bowers fc Deeming (1984), 


Islam et al.l (1201511 calculated the allowed range 
of orbital inclinations of the sys tem. Assuming 


a sta r with spectral type 09 I (IBodaghee et al 


201211 . t he orbital inclinat ion was constrained to 
39-57° ( Islam et al.l 2015 1. A main-sequence B- 


type star yields orbital inclinations between 60- 


77° (jBodaghee et al.ll2012i llslam et al.N2ni5[l . 


We deriye a 4.2378±0.0004d orbital period for 
IGR J16393-4643 using a DFT, which is co n sisten t 
with the results from iGorbet fc Krimnil ( 201,31 1. 
Using an O — C analysis (see Section 12.21) . this 
is further refined to 4.23810±0.00007d. We note 
that we obtain a bad fit in the mid-eclipse time 
between MJD 56079-56745. As a result, we only 
use data between MJD 53416-56078 in our O — C 
analysis (see Figures [THSl- [3]). Using the quadratic 
orbital change function (see Equation 2]) , we find 
the orbital period deriyatiye to be -5±4xl0“^d 
d“^, which is consistent with zero. The dura¬ 
tion of the observed eclipse was calculated to be 
31 I 5 ks (0.36tao6 yielding an eclipse half-angle 
of 15l3° (see Table [I]). We find these t o be c onsis- 
tent with the result from llslam et al.l (|2015^ . The 
source flux does not reach 0 counts cm“^ s^ in the 
folded light curves during eclipse (see Figure |4]). 
We interpret this dip as an eclipse since the fea¬ 
ture is persistent over many years of data. The 
rapid ingress and egress requires obscuration by 
clearly defined boundaries that are suggestive of 
an o bject such as the mass donor in the system 
('e.g. IColev et al.ir2014ll . We discuss the nature of 
the non-zero flux during eclipse in Section 14.21 


We calculate the predicted eclipse half-angle ©e 
as a function of inclination angle of the system (see 
Figure [S]). The calculation assumes a neutron star 
mass of 1.4 Mq, and the primary stellar masses 
and radii given in Table [6l We calculate the min¬ 
imum inclination angle of the system, Zmin, that 
is consistent with the measured eclipse half-angle 
(see Table HI). We find that stars with spectral 
types BO V, BO-5 III and BO I satisfy the con¬ 
straint imposed by the eclipse half-angle (see Ta¬ 
ble [6]) . We note that while a B5 III star satisfies 
the constraint imposed by the minimum value of 
the eclipse half-angle under the assumption that 
the neutron star is 1.4 Mq, this spectral type does 
not satisfy the eclipse half-angle for a more mas¬ 
sive neutron star. 


Since the X-ray luminosity is lower than 
what would be exp ected for Roche-lobe overflow 
( Kaoer et al.l 12004!) . we can attach an additional 


constraint assuming that the mass donor underfills 
the Roche-lobe radius (see FigureH]). We note that 
while the derive d masses and radii for th e spectral 
type BO I from ICarroll fc Ostlid (I 2 OO 6 I) satisfies 
the eclipse half-angle, the assumed rad ius would 
be larger than the Roche-lobe radius (jEggleton 
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Table 6 

Physical Parameters for Previously Proposed Mass Donors for IGR J16393-4643 


Spectral Type 

M/Mq 

qa 

R/Rq 

Rl/R0^ 


BO III 

20 

0.070 

13 

18.5 

68 

B5 III 

7 

0.200 

6.3 

11.7 

79 

BO V 

17.5 

0.080 

8.4 

17.5 

79 

BO I 

25 

0.056 

25'^ 

20.4^ 

41 

Note.—T he 

values 

in italics are 

obtained 

from 


“ The mass ratio, g, is defined as M^/Mc where is the com¬ 
pact object and Me is the donor star. 

^ The definition for the Roche lobe, Rl, as given in lEggleton 
( 1983h . assuming Mns is 1.4 M 0 . 


The minimum inclination angle of the system that is consistent 
with the measured eclipse half-angle. 

A BO I classification significantly overfills the Roche-lobe and 
are therefore, it is excluded from our analysis. 


19831 1. Therefore, an BO I spectral type must be 
excluded (see Table HI). 

3.2. IGR J16418-4532 

IGR J16418-4532 is a candidate Supergiant 
Fast X-ray Transient (SFXT ) first discovered with 
the INTEGRAL satellite bv lTomsick et ahl ( 2004h 


at a flux of 3x10“^^ erg cm“^ s“^ in the 20-40keV 
band. The near-infrared spectral energy distribu¬ 
tion of the most probable Two Micron All Sky 
Survey (2MASS) counterpart was measured with 
the 3.5 m New Tech nology Telescope (N TT) at 
La Silla Observatory, iRahoui et al.1 (l2008ll fo und a 
spectral type of 08.5 I. Coleiro et al.l (l201^ pro¬ 
posed a spectral type of BN0.5 la based on features 
in the near-infrared spectrum such as Br(7 — 4) 
and the emission and absorption of neutra l helium. 
Using XMM-Newton, Walter et al. ( 2006[l found a 
1246±100s modulation, which was interpreted as 
a neutron star rota tion period. This was later re¬ 


fined to 1212±6s (ISidoli et al.l 1201211. also using 


XMM-Newton. Recently, iDrave et al.1 (|2013^ fur¬ 
ther refined the rotation period to 1209.1±0.4s 
with XMM-Newton. A ~3.73d orbital period 
was found using data from the Swift BAT and 
the RXTE ASM instruments, where Porb was 


reported as 3.753±0.004d and 3.7389±0.0004d, 

2n0bll. 


respectively (ICorbet et al.1 


tended ASM dataset, iLevine et al.l ( 2011 1 found 


Using an ex- 


a 3.73886±0.00028d period , which is co n sisten t 
with the earlier r e sult fr om ICorbet et al.l ( 20061 1 . 
Corbet fc KrimrnI ( 2013ll further refined this to 
3.73886±0.00014d using BAT. A -14.7d mod¬ 
ulation was found using BAT and INTEGRAL, 
which was interpreted as a superorbital period 
( Corbet fc Krimm 20131 : IPrave et al. 2013 1. 

The ASM and BAT light curves folded on 
the orbital period revealed a sharp dip with 
near zero mea n flux, which was interpreted as 


a total eclipse (ICorbet et al. 
observations of the e clipse 


ray Telescope (XRT) (iRomano et al 


2006ll . Subsequent 
included Swif t X- 
2012 1 and 


INT EGRAL (iDrave et al l l2013ll . With Swift 
BAT, Romano et al.l ( 2012h constrained the eclipse 
half-angle to be 0.17±0.05 of the orbital pe¬ 
riod, corresponding to a duration of 0 . 6 ± 0 . 2 d 
(55±16ks). The duration of the eclipse was 
found to be ~0.75d in the archival data set 
from INTEGRAL/IBI S, covering t he time pe¬ 
riod MJD 52650-55469j|D^^re^[20ll). While 
the estimate in IDrave et al.l (l2013li ' is significantly 


l arger than the constraints from iRomano et al 
(I2OI2I I. a lower limit of ^ 0 . 583 d was found ir 
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Fig. 5.— The black curves show the predicted 
eclipse half angle of IGR J16393-4643 as a function 
of inclination angle for stars with the indicated 
spectral types. The red and black dashed lines 
indicate the eclipse half angle and estimated er¬ 
ror as measured by Swift BAT. We assume a neu¬ 
tron star mass of 1.4 Mq (top) and of mass 1.9 Mq 
(bottom) and typical masses and radii for the as¬ 
sumed companion spectral type (see Tabled). The 
blue vertical dashed lines indicate the lower limit 
of the inclination angle. Inclinations to the left of 
these correspond to stars that overfill the Roche- 
lobe. 


a combined st udy with IN TEGRAL/IBIS and 
XMM-Newton ( Drave et al.ll201^ . 

We derive a 3.73863±0.00015d orbital pe¬ 
riod for IGR J16418-4532 using the fundamen¬ 
tal peak in the power spectrum, while the first 
harmonic yields 3.73882±0.00011d. Using an 
O — C analysis (see Figures [BH), we refine this 
to 3.73881±0.00002d. Using the quadratic orbital 
change function (see Equation SI), we find the 
orbital period derivative to be 0.7±1.0xl0“^d 
d“^, which is consistent with zero. Folding the 
light curve on our refined orbital period (see Fig¬ 
ure III), we calculate the duration of the observed 
eclipse to be 57±Iks (0.66±0.0Id). This yields 
an eclipse half-angle of 3I.7lo;g° (see Table [2]). 
We note that we hnd the eclipse half-angle to be 
3I.5±0.6° assuming a symmetric step-and-ramp 
function. We find the eclips e properties to be con ¬ 
siste nt with the r e sults f rom iRomano et ah ( 2012 1 


and iDrave et al.l (|2013[ l. Under the assumption 


that the mass and radius of the proposed mass 


Fig. 6.— Log-log plot of stellar mass as a func¬ 
tion of stellar radius for IGR JI6393-4643. The 
dark shaded region indicates the allowed spec¬ 
tral types that satisfy the constraints imposed by 
both the eclipse duration and Roche-lobe size. 
The light shaded region only indicates spectral 
types that satisfy the observed eclipse duration. 
Stellar masses and radii are reported in Table [SI 
The red, green and blue lines indicate interpola¬ 
tions for main-sequence, giant and supergiant lu¬ 
minosity classes, respectively. The dashed, dot¬ 
ted and dash-dotted lin es indicate stellar masse s 
and radii der i ved fr om Carroll fc Ostlid (2006), 


[Martins et al. (2005) and Allen ( 2000l l. 


respec¬ 


tively. 


dono r are 31.54 Mq and 21.41 Rq ([Martins et al 


20051) . the duration of the observed eclipse is con¬ 


sistent with the proposed mass donor, where we 
find the orbital inclination to be 60-63° for an 
08.5 I spectral type (see Figures [H and [5|). 

Since the X-ray luminosity is lower than 
what would be exp ected for Roche-lobe overflow 
( Kaoer et al.l 120041) . we can attach an additional 


constraint assuming that the mass donor under¬ 
fills the Roche-lobe radius (see Figure[n|). We note 
that while the deri ved masses and radii for the 
spectral types from [Martins et al.l ( 2005ll satisfy 
the eclipse half-angle, the assumed rad ius would 
be la rger than the Roche-lobe radius ( EggletoJ 


I983I) . Therefore, an 08.5 I spectral type must be 


excluded. 
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Table 7 

Physical Parameters for Previously Proposed Mass Donors for IGR J16418-4532 


Spectral Type 

M/Mq 

qa 

R/Rq 

Rl/R0^ 

i° 

08.51 

31.54 

0.044 

21.4T 

20.71"= 

61-63 

08.51 

33.90 

0.041 

22.2(F 

21.35"= 

60-62 

091 

29.63 

0.047 

21.76’^ 

20.17"= 

58-60 

091 

31.95 

0.044 

22.6(f 

20.82"= 

57-59 


Note. —The values in italics are obtained fromlMartins et al.l 


(120051) . 


“ The mass ratio, q, is defined as M^/Mc where is the com¬ 
pact object and Me is the donor star. 

^ The definition for the Roche lobe, i?L, as given in EggletonI 
(llOS.lh . assuming Mns is 1.4 M 0 . 

^ These spectral types significantly overfill the Roche-lobe and 
are therefore excluded from our analysis. 


3.3. IGR J16479-4514 

IGR J16479-4514 is an intermediate SFXT, 


(Chatv et al. 

2008 

R.ahoui et al. 2008h 

lab (Nesooli et al. 

200§l) mass donor. 


First dis- 

cover ed by the INTEGRA L satellite in 2003 Au¬ 


gust (IMolkov et al.l 1200311 . the fluxes in the 18- 


25keV and 25-50 keV energy bands were found 
to be ~12mCrab and ~8mCrab, respectively. 
The flux was later shown t o increase by a fac - 
tor of ~2 on 2003 Augus t 10 (iMolkov et al.ll2003ll . 


Using Swift BAT data, I Jain et ah ( 2009 ) found 


the presence of a 3.319±0.001 day modulation, 
which was interpreted as the orbital period. A 
3.3193±0. 0005dav modu l ation was independently 
found by Romano et al. ( 20091) also using BAT. 
Corbet fc Kriminr ( 2013 ) found the orbital period 
to be 3.3199 ±0.0005d, which i s con sistent with the 
result s from Jain et al. ( 20091 ) and Romano et al. 
( 20091 ). The presence of a 11.880±0.002d su¬ 


perorbital period was found by ICorbet fc Krimm 
( 2013 ) using BAT. Drave et al. ( 2013h reported a 
11.891±0.002d superorbital period using INTE¬ 
GRAL, confirming the result. No pulse period has 
been identified. 

The BAT light curves folded on the orbital 
period revealed a sharp dip, which was inter¬ 
preted as an eclipse with a proposed duration of 


^52 ks ( Jain et al. 20091) . This confirmed an ear¬ 
lier XMM-Newton observation where a decay from 
a higher to lower flux st ate was interprete d as 
the ingress of an eclipse ( Bozzo et al. I I 2 OO 9 I) . A 
2012 Suzaku observation covered ^80 percent of 
the orbital cycle of IGR J16479-4514, where tem¬ 
poral and spectral properties were analyzed dur - 
ing eclipse and out-of-eclipse (|Sidoli et al.ll2013l) . 
Since the ingress of the eclipse was not covered 
in Suzaku observation, the exact duration of the 
eclipse could only be constra ined to 46-143 ks 
(0.53-1.66d) (ISidoli et al.ll2013l) . 

Using a DFT, we derive the orbital period of 
IGR J16479-4514 to be 3.31998±0.00014d. We re¬ 
fine this to 3.3I96I±0.00004d using an 0 — C anal¬ 
ysis (see Figures [TH3|) and fold the light curve on 
our refined orbital period to calculate the eclipse 
half-angle (see Figure [ini). Using the quadratic 
orbital change function (see Equation |4]) , we find 
the orbital period derivative to be 3±2xl0“^d 
d“^, which is consistent with zero. We calcu¬ 
late the duration of the observed eclipse to be 
52±3ks (0.60±0.03d) , whic h is consistent with re¬ 
sults from ljain et al.l ( 20091) . This yields an eclipse 
half- angle of 31. 9^^ ' ^° (se e Tablellt. Using values 


( 20051 ) for the masses and radii 
of the proposed spectral type of the mass donor, 
the duration of the observed eclipse is consistent 
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Fig. 7.— Swift-BAT light curve of IGR J16418- 
4532 in the 15-50keV band folded on the or¬ 
bital period (top) using 20 bins. TO is defined 
at BMJD 55087.721, corresponding to mid-eclipse. 
A detailed folded light curve with 80 bins (bot¬ 
tom) is fit with both a symmetric “step and ramp” 
function (green) and asymmetric “step and ramp” 
function (red), which model the eclipse. The sym¬ 
metric “step and ramp” function was shifted ac¬ 
cordingly. 


with the proposed mass donor (see Figures ITT]) , 
where we hnd the orbital inclination to be 54- 
58° or 47-51° for a 08.5 1 or 09.5 lab spectral 
type, respectively (see Table B- While the previ¬ 
ously preposed spectral types satisfy the eclipse 
half-angle, we note that the radius of the pro- 
posed spectral t ype is larger than the Roche-lobe 
( Eggletonlll983 ). Therefore, the previously pro¬ 
posed 08.5 1 or 09.5 lab spectral types mu st be 
excluded (see Figure IT^. ISidoli et al.l ( 201 ^ pro¬ 
posed a mass donor with mass ^35 Mq and ra¬ 
dius ~ 2 Oi? 0 , which could satisfy the constraints 
imposed by the Roche-lobe radius. 


3.4. IGR J18027-2016 (=SAX J1802.7- 
2017) 

IGR J18027-2016 (=SAX J1802.7-2017) is 

josed to host 

imnli 


either 

a B1 Ib (i 

Urreion et al. 

B1 I ( 

Mason et al. 

20 Ilf) mass 


or BO- 
donor. First 


detected with Bep poSAX in 2001 September 
( Augello et al. 2003ll . the average flux in the 0.1- 
10 keV energy band was found tobe3.6x10“^^ ergs 
cm“^ s“^. Pulse-timing analysis suggested a 


Fig. 8 .— The black curves show the predicted 
eclipse half angle of IGR J16418-4532 as a function 
of inclination angle for stars with the indicated 
spectral types. The red and black dashed lines 
indicate the eclipse half angle and estimated er¬ 
ror as measured by Swift BAT. We assume a neu¬ 
tron star mass of 1.4 Mq (top) and of mass 1.9 Mq 
( bottom) and typical masses and radii for the as¬ 
sumed companion spectral type (see Table[7]). The 
blue vertical dashed lines indicate the lower limit 
of the inclination angle. Inclinations to the left of 
these correspond to stars that overfill the Roche- 
lobe. 


~4.6day orbital period ( Augello et al.l l2003ll . 
which was later refined to 4.5696±0.0009days 
INTEGRAL (iHill et al.l l2005ll. Gombin- 


usmg 


i ng th e mid -eclipse times der ived in I Augello et al.l 


([20031) and iHill et al.l (j2005r i with later observa¬ 
tions with Swift BAT and INTEGRAL in an 
O — C analysis, the orbit al period was fu rther 
refined to 4.5693±0.0004d ( Jain et al. 2009ll . Fit¬ 
ting a quadratic model to the derived mid-eclipse 
times, a per iod derivative of (3.9±1.2) x 10“^ d d“^ 
was found (l.lain et al.ll27ioil . Using INTEGRAL, 
Falanga et al.l ( 2015[ l recently refined the orbital 


period and period derivative to 4.5697 ±0.0001d 
and (2 .1±3.6)xl0~^dd~^, respectively. lAugello et al 
( 2 OO 3 I I found the presence of a ^139.6 s modula¬ 
tion using BeppoSAX, which t hey interpeted as 


the neutron star rotation period. iHill et al.l (12005^ 


found the pulse period using XMM-Newtoii to be 
139.61±0.04 s, which is co n sisten t with an earlier 


result from lAugello et al.l (120031 1. However, no 


evidence for superorbital modulation was found 
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Table 8 

Physical Parameters for Previously Proposed Mass Donors for IGR J16479-4514 


Spectral Type 

M/Mq 

qa 

R/Rq 

Rl/R©'’ 

i° 

08.51 

31.54 

0.044 

21.4T 

19.13'= 

55-58 

08.51 

33.90 

0.041 

22.2(T 

19.73" 

54-56 

09.51 

27.83 

0.050 

22.1T 

18.14" 

48-51 

09.51 

30.41 

0.046 

23. IT 

18.84" 

47-49 


Note. —The values in italics are obtained fromlMartins et al.l 


(120051) . 


“ The mass ratio, q, is defined as M^/Mc where is the com¬ 
pact object and Me is the donor star. 

^ The definition for the Roche lobe, i?L, as given in EggletonI 
(Il08,lh . assuming Mns is 1.4 M 0 . 

^ These spectral types significantly overfill the Roche-lobe and 
are therefore excluded from our analysis. 


(jCorbet fc Krimmll2013r) . 


Augello et al.l (120031) derived the epoch of the 


NS superior conjuction to occur during the first 
50 ks of an observation where IGR J18027-2016 
was undetected. Since the epoch of superior 
conjuction corresponds to a time where the mid- 
eclipse is expected to occur, Augello et al. ( 2003h 


suggested the presence of a possi ble eclipse with a 


half- duration of ( ).47± 0.10davs. iHill et al 
and I Jain et ahl ( 20091) later confirmed 


(1200, 5h 
this re¬ 
sult where the eclispe half-angle was found to be 
0.61±0.08radians (34.9il i4.6 °) and '^0. 604radians 
(~34.6°), respectively. IPalanga et al. ( 2015 ) re¬ 
cently refined the eclipse half-angle to 31 ±2 ° using 
INTEGRAL. 

Pulse-timing and radial velocity curves have 
helped place constraints on the physical proper¬ 
ties on both the donor star as well as the compact 
object. Using a pulse-timing analysis with Bep- 
poSAX and XMM-Newton, the projected semi¬ 
major axi s of the neutron star was found to be 
68 ±llt-s ( Hill et ^ 2005 ). Assuming a 1.4 Mq 
neutron star. Hill et aH ( 20051) constrained the 


mass and radius of the mass donor star to 18.8- 
29.3 Mq and 14.7-23.4 Rq, re spectively. Usin; 


Swift BAT and INTEGRAL, Jain et al. (200' 


further constrained the radius of the donor star 
to 16.4-24.7i?0. The mass donor was observed 


between 2010 May 26 and 2010 September 8 
(MJD 55342-55447) with the European South¬ 
ern Observatory Very Large Telescope {ESO 
VLT, iMason et ahl 1201 ll) . The semi- amplit ude 
of the radial velocity of the mass donor, Kg 


was found to be 23.8±3.1km s ^ (^Mason et al 


20111) . Since the projected semi-major axis of 


the neutron star can be expressed in terms of a 
radial velocity semi-amplit ude, ATx, tbe ratio be¬ 
tween the masses of the compact object and the 
mass don or can be calcula ted according to Equa - 


tion 2 in Ash et al. ( 19991) . 


Mason et al.l ( 201ll) 
found the mass ratio, q, to be 0.07±0.01. Using 
the mas s ratio and the eclipse half-angle mea- 
( 1200 , 5h . the mass and radius 


sured in Hill et al 


of the donor star were refined to values between 
18.6±O.8M0 and 16.8±1.5i?0 at edge-on incli¬ 
nations to 21.8±2.4M0 and 19.8± 0.7R(7^ whe re 


the donor star fills the Roche-lobe (Mason et al 


20111 ). The mass of the compact object was con¬ 


strained to be between 1.36±0.21 and 1.58±0.27 Mq 
in the two limits. The large error on the esti mate 
of the eclipse half-angle from iHill et al. ( 2005 ) con¬ 
tributes significantly to the uncertainties on these 
measureme nts. Using a s i milar analysis with IN¬ 
TEGRAL, iFalanga et al. ( 201,51 ) constrained the 
mass of the neutron star to 1 . 6 ±O. 3 M 0 . 

We derive a 4.57022±0.00013d orbital period 
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Table 9 

System Parameters for IGR J18027-2016 


Parameter 

Value 

Rorb“ 

4.56993±0.00003d 

Hpulse'' 

139.61±0.04s 

flx smi** 

68 ±llt-s 

Ro" 

23.8±3.1km s”! 

Tjnid 

BMJD55083.82±0.01 

f{Mf 

16±1 Mq 


0.07±0.01 

0e 

34±2° 

-^^donor 

18.6±0.9-19.4±0.9Mo 

-^donor 

17.4±O.9-19.5tEJ-?i?0 

a 

31.4-33.2 Rq 

Mns 

1.37±0.19-1.43±0.20Mq 

i 

73.3-90° 


Note. —“ The orbital period is refined 
using an O — C analysis. 

^ The pulse period, projected semi-major 
axis and mass function are given in 


Hill et al.l (|2nn,^ . 


° The semi-amplitude of the radial veloc¬ 
ity of th e mass donor and m ass ratio are 
given in 


Mason et al.l (120 llh . 
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Fig. 9.— Log-log plot of stellar mass as a func¬ 
tion of stellar radius for IGR J16418-4532. The 
dark shaded region indicates the allowed spectral 
types that satisfy the constraints imposed by both 
the eclipse duration and Roche-lobe size. The 
light shaded region only indicates spectral types 
that satisfy the observed eclipse duration. Stel¬ 
lar masses and radii are reported in Table [71 The 
dashed, dotted and dash-dotted li nes indicate stel- 
lar masses and radii derived from Carroll fc Ostlie 


( 2006il . [Martins et al. ( 2005 1 and Allen! ( ^OOll . re- 
spectively. 


for IGR J18027-2016 usin g a DFT, which is consis- 


g a J 

tent with the results from jHill et al.l ( 2005h . Using 
an O — C analysis (see Figures (THS]), we refine this 
to 4.56993±0.00003d. Using the quadratic orbital 
change function (see Equation |T]) , we find the or¬ 
bital period derivative to be 0.2±l.lxl0“^d d“^, 
which is consistent with zero. Folding the light 
curve on our refined orbital period(see Figure [T51) . 
we measure the duration of the observed eclipse 
to be 74±4ks (0.85±0.05d). We find the eclipse 
half-angle to be 34±2° (see Table |T]). Since the 
X-ray luminosity of IGR J18027-2016 is modest, 
we again attach the constraint that the donor star 
underfills the Roche-lobe (see Figures [T4l and (TSl) . 
This constrains the mass and radius of the donor 
star as well as the mass of the neutron star (see 
Figure dsi). Using the eclipse half-angle and the 
expression for th e mass donor when the mass ratio 
is known (Eq.4; [Mason et al.l 201 ih , we calculate 
the donor star mass and radius as well as mass 
of the compact object. We find the mass and 
radius of the donor star to be 18.6±O.9M0 and 


7 0.0015 
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Fig. 10.— Swift-BAT light curve of IGR J16479- 
4514 in the 15-50 keV band folded on the or¬ 
bital period (top) using 20 bins. TO is defined 
at BMJD 55081.57, corresponding to mid-eclipse. 
A detailed folded folded light curve with 80 bins 
(bottom) is fit with both a symmetric “step and 
ramp” function (green) and asymmetric “step and 
ramp” function (red), which models the eclipse. 
The symmetric “step and ramp” function was 
shifted accordingly. 

17.4±O.9i?0 and 19.4±O.9M0 and lO.Sto^yR© 
the two limits (see Table [3]). In the allowed limits, 
we constrain the mass of the neutron star to be¬ 
tween 1.37±O.19M0 and 1.43±O.2OM0 (see Ta¬ 
ble in|) • While our results a r e in a greement with 
calculations in Mason et al.[ ( 201l!l . we note that 
the error estimate is only marginally improved in 
our analysis. The driving factor on the error esti¬ 
mate of the neutron star mass is the uncertainty 
of ~13% on the semi-amplitude of the radial ve¬ 


locity of the donor star as reported in iMason et al 

( 1201111 . 
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Fig. 11.— The black curves show the predicted 
eclipse half angle of IGR J16479-4514 as a function 
of inclination angle for stars with the indicated 
spectral types. The red and black dashed lines 
indicate the eclipse half angle and estimated er¬ 
ror as measured by Swift BAT. We assume a neu¬ 
tron star mass of 1.4 Mq (top) and of mass 1.9 Mq 
(bottom) and typical masses and radii for the as¬ 
sumed companion spectral type (see TablejS]). The 
blue vertical dashed lines indicate the lower limit 
of the inclination angle. Inclinations to the left of 
these correspond to stars that overfill the Roche- 
lobe. 


Fig. 12.— Log-log plot of stellar masses as a 
function of stellar radius for IGR J16479-4514. 
The light shaded region indicates the allowed spec¬ 
tral types that satisfy the observed eclipse. The 
dark shaded region indicates the spectral types 
where both the eclipse duration and Roche-lobe 
radius constraints are satisfied. Stellar masses and 
radii are reported in Table [5] The dashed, dot¬ 
ted and dash-dotted lin es indicate stellar masse s 
and radii der i ved fr om Carroll fc Ostliel ( 2006 1. 
iMartins et al. (2005) and Allen (2000), respec¬ 
tively. 


3.5. XTE J1855-026 

XTE J1855-026 is an SGXB disco vered during 
RXT E scans along the Galactic plane (jCorbet et al 
I 999 I I. Through 11 scanning observations of the 
Scutum arm, the 2-10 keV flux of XTE J1855- 
026 varied from an u pper limit of 10 cou nts s“^ 
to 136±15counts s“^ ( Corbet et al.lll999l) . Using 
the RXTE P r oport ional Counter Array (PGA), 


Corbet et al. ( 1999ll found a 361.1±0.4s mod¬ 


ulation, which was int erpreted as the n eutron 


star rotation period. ICorbet fc Mukail (120021) 


later refined the pulse period to 360.741±0.002s 
with the PGA. An analysis using the RXTE 
All-sky Monitor (ASM) revealed the presence of 
a 6.0724±0.0009d modulation, which is inter¬ 


prete d as the orbital period ([Corbet fc Mukai 
2002). Using an O — C analysis with INTE¬ 
GRAL, Falanga et al. ( 20151) recently refined this 
to 6.07415±0.00008d. No significant orbital pe¬ 
riod derivative was fo und using the quadr atic or¬ 
bital change function ( Falanga et al.l 12015 ). This 


places XTE J1855-026 in the wi nd-fed super ¬ 
giant region of the Corb et diagram ( Corbet 19861) . 


Corbet &: Mukail (|2002l) constrained the eccentric¬ 


ity of XTE J1855-026 to e<0.04 using pulse-timing 
analysis. The projected semi-major axis of the 
neutron star was found to be 82.8±0.81t-s for 
a circular orbit and 80.5±1.41t-s for the eccen¬ 


tric solution ([Corbet fc Mukail l2002l) . In this sec¬ 
tion, we consider the scenario where the orbit is 
circular. We calculate constraints on the mass 
and radius of the donor star in the more com¬ 
plicated scenario where the orbit is considered to 
have a modest eccentricity in Section 14.31 Using 
optical and near-infrared spectroscopy obtained 
with the William Her schel Telescope (WHT), 
Negueruela et al. ( 20081 ) found the mass donor 


to be a supergiant with spectral type BO laep. 

The light curves folded on the orbital period 
reveal a sharp dip, which was interpreted as an 
eclipse with a total pha se duration of 0.198-0.262 
( Corbet fc Mukail 2002 ). The phase duration cor¬ 
responded to an eclipse half-angle 36°<0e<47°. 
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Fig. 13.— Swift-BAT light curve of IGR J18027- 
2016 in the 15-50 keV band folded on the or¬ 
bital period (top) using 20 bins. TO is dehned 
at BMJD 55083.82±0.01, corresponding to mid- 
eclipse. A detailed folded light curve with 80 bins 
(bottom) is ht with both a symmetric “step and 
ramp” function (green) and asymmetric “step and 
ramp” function (red), which model the eclipse. 
The symmetric “step and ramp” function was 
shifted accordingly. 


The eclipse duration was found to be 93±3ks 
(1.08 ±0.03d) in th e archival INTEGRAL data 
set ( Falanga et al.l 2015[l . corresponding to an 
eclipse half-angle of 32±1°. Thi s measurement is 
somew hat lower than result from ICorbet fc Mukai 

(|2002ll . 

Optical radial velocity curves recently obtained 
with the Isaac Newton Telescope (INT), the Liver¬ 
pool Telescope (LT) and the WHT help place ad- 
ditio nal constraints on the components of the sys¬ 
tem ( Gonzalez-Galanll2015 b The semi-amplitude 
of the radial velocity of the donor star was found 
to be 26.8±8.2km s“^. Expressing the projected 
semi-major axis of the n eutron star as a radia l 
velocity semi-amplitude, Gonzalez-GalanI ( 2015h 
found the ratio between the masses of the com¬ 
ponents of the system to be 0.09±0.03 and notes 
that a large eccentricity of ^0.4-0.5 was found in 
the optical orbital solutions. This strongly con- 
trasts with the X - ray o rbital solution reported in 


Gorbet fc Mukail (|2002[ 1. suggesting that caution 


must be t aken in interpre t ing th e optical orbital 
solutions. iGonzalez-GalanI (12015!) rehned the spec¬ 
tral type of the mass donor to a BN0.2 la super- 


Fig. 14.— The black curves show the predicted 
eclipse half angle of IGR J18027-2016 as a function 
of inclination angle for stars with the indicated 
spectral types. The red and black dashed lines 
indicate the eclipse half angle and estimated error 
as measured by Swift BAT. 


giant and found the mass and radius of the donor 
star to be Mq and 27tiioRe, respectively. 

We derive a 6.07411±0.00014d orbital period 
for XTE J1855-026 using the fundamental peak in 
the power spectrum. Using an O — C analysis (see 
Figures we rehne this to 6.07413±0.00004d 
and fold the light curve on the rehned orbital 
period to calculate the eclipse half-angle (see Fig¬ 
ures m and ini). We calculate the duration of 
the observed eclipse to be 98±2ks (1.13±0.02d), 
yielding an eclipse half-angle of 33.6±0.7degrees 
(see Table [ini). This indicates that our derived 
eclipse du ration is somewh a t less than the re¬ 
sult from iGorbet fc Mukai ( 2002 ) and is con- 


sistent with the measurement in iFalanea et al 


( 20151 ). Using the quadratic orbital change func¬ 
tion (see Equation SI) , we hnd the orbital period 
derivative to be 0.0±0.5xl0“^d d“^, which is 
consistent with zero. 

Since the upper limits of the stellar mass and 
radius are constrained by the orbital inclination 
where the Roche lobe is just hhed, we again at¬ 
tach constraints on the mass and radius of the 
donor star as well as the mass of the compact 
object. We hnd the inclination where the donor 
star hlls the Roche lobe to be 76.4°. The stellar 
mass and radius of the donor star are found to be 
19.6±1.1M0 and 21.5±0.5Ro and 20.2±1.2Mo 
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Fig. 15.— Log-log plots of stellar masses as a 
function of stellar radii for IGR J18027-2016. The 
shaded region indicates the allowed spectral types 
to satisfy the observed eclipse duration and pulse¬ 
timing constraints. Stellar masses and radii are 
reported in Table [9l Stars and crosses indicate 
the sp ectr al types according t o ICarroll fc: Ostlie 


me sp ectr a .1 types according t o it^arroii U 
( 200611 and Lefever et al. ( 2007ll . respectively. 


Fig. 16.— Swift-BAT light curve of XTE J1855- 
026 in the 15-50 keV band folded on the orbital 
period (top) using 20 bins. TO is defined as 
MJD 55079.0685, corresponding to mid-eclipse. A 
detailed folded folded light curve with 80 bins 
(bottom) is fit with both a symmetric “step and 
ramp” function (blue) and asymmetric “step and 
ramp” function (red), which models the eclipse. 


and 23.O±O.5i?0 in the two limits (see Table [T51) . 
We hnd the mass of the neutron star to be between 
1 .77TO.55M0 and 1.82±O.57M0 (see Table (HI), 
where the driving factor on the large error esti¬ 
mate is the uncertainty of ^31% on the radial ve- 
lo city semi-amplitude of t he donor star as reported 


Gonzalez-GalanI (|2015[ 1. 


4. Discussion 

We discuss our findings of IGR J16393-4643, 
IGR J16418-4532, IGR J16479-4514, IGR JI8027- 
2016 and XTE J1855-055. The radii for the previ¬ 
ously proposed spectral types in IGR JI64I8-4532 
and IGR J16479-4514 would significantly overfill 
the Roche-lobe, which suggests an earlier spectral 
type. Below, we discuss in detail the nature of 
the mass donors in each system, mechanisms that 
could result in the residual emission observed in 
IGR J16393-4643 and comment on the nature of 
the eclipse profiles. 


4.1. What is the nature of the mass donors 
in each system? 

11.1. IGR J16393-4643 

Our results show that stars with spectral type 
BO V and BO-5 III satisfy the constraints imposed 
by the eclipse half-angle as well as the Roche-lobe 
(see Section EH). While some supergiant stars 
such as a BO I satisfy the eclipse half-angle, the re¬ 
quired radius would be larger than the Roche-lobe. 
We calculated the Roche-lobe radius for stars of 
spectral type BO I to be 18.5 Rg^ . This is clearly 
smalle r tha n the radii rep orted in ICarroll fc Ostlie 
( 20061 ) and lAllenI ( 200(il) . which are 25 Rq and 
30 Rq, respectively (see Table E]). Since the ra¬ 
dius for this proposed supergiant is too large to 
satisfy the constraint imposed by the Roche-lobe, 
we suggest that if the donor star is a supergiant 
then it must be a slightly earlier spectral type. 

IGR J16393-4643 was observed in the near- 
infrared on 2004 July 9 (MJD 53195.3) using 
the 3.5 m New Techn ology Telescope (N TT) at 
La Silla Observato ry ( Chatv et al.l 12001^ . While 
Ghatv et al.l ( 20081 ) conclude that the spectral type 
of the donor star is BIV-V based on the spectral 
features and spectral energy distribution (SED), 
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Table 10 

System Parameters for XTE J1855-026 


Parameter 

Value 

^’orb“ 

6.07413±0.00004d 

impulse'’ 

360.741±0.002s 

-Ppulse'’ 

1.5±3.6s s-^xlO"® 

Ox sini^ 

82.8±0.81t-s 


26.8±8.2km s-^ 

Tjnid 

MJD55079.07±0.01 

f{Mf 

16.5±0.5Mq 

q- 

0.09±0.03 

0e 

33.6±0.7° 

-^^donor 

19.6±1.1-20.2±1.2Mq 

-^donor 

21.5±O.5-23.O±O.5i?0 

a 

38.9-39.2 i^Q 

Mns 

1.77±O.55M0-1.82±O.57M0 

% 

76.4-90° 


Note. —“ The orbital period is refined us¬ 
ing an O — C analysis. 

^ The pulse period, derivative of the pulse 
period, projected se mi-major axis and mas s 
function are given in Corbet fc Mukail ( 20021) . 

The semi-amplitude of the radial velocity of 
the mass donor and m ass ratio are given in 
Gonzalez-GalanI ( 2015h . 


21 





















Fig. 17.— The black curves show the predicted 
eclipse half angle of XTE J1855-026 as a function 
of inclination angle for stars with the indicated 
spectral types. The red and black dashed lines 
indicate the eclipse half angle and estimated error 
as measured by Swift BAT. 


we note that Nespoli et al. ( 2008ll proposed the 
donor star to be a K or M supergiant using the 
same SED. Eurthermore, observations with the 
Chandra observatory show that the previously 
proposed counterpart is positionally incons istent 


with the X-ray source (IBodaghee et al.ll2012 ). Us- 
ing the GLIMPSE survey, Bodaghee et al. ( 2012 1 


proposed that the counterpart must be a distant 
reddened B-type star. 

4.f.3. IGR JI 64 I 8-4532 

Our results show that while the derived masses 
and radii for the previously proposed spectral 
types satisfy the eclipse half-angle (see Table III), 
the radius would b e larger than the Roche-lobe 
size ( EggletonlllQS.lI l An 08 I star has a ma ss of 
28 Mp acco rding to ICarroll fc Ostlie ( 2006ll and 
AllenI ( 20001) . We find the maximum radius for 


a 28 Mq star to satisfy the constraint imposed 
by the Roche-lobe to be 18.2 Rf7> . This is clearly 
smalle r tha n the radii rep orted in Carroll fc Ostlid 
( 2006ll and lAllenl ( 2000 ll . which are 22 Rq and 
2 Oi?0, respectively (see Table [T]). Since the ra¬ 
dius for each proposed spectral type is too large 
to satisfy the constraint imposed by the Roche- 
lobe, it is our contention that the donor must be 
an earlier spectral type. 

We find that spectral classes 07.5 I and ear- 


Fig. 18.— Log-log plots of stellar masses as 
a function of stellar radii for XTE J1855-026. 
The shaded region indicates the allowed spec¬ 
tral types to satisfy the observed eclipse dura¬ 
tion, pulse-timing and optical radial velocity con¬ 
straints. Stellar masses and radii are reported in 
Table [TUI Stars and crosses indicate t he sp ectral 
types according to ICarroll &: Ostlid ( 2006ll and 


Lefever et al. ( 2007ll . respectively. 


lier satisfy the constraint imposed by the Roche- 
lobe radius (see Table fTTl) . The ratio between the 
radius of the donor star and that of the Roche- 
lobe, /3, is found to exceed 0.9, which is con- 
sistent with other HMXB s that host supergiants 
( Joss fc R,appaport 1984l l. Transitional Roche- 


lobe accretion has been proposed in IGR J16418- 
4532 where a fraction of the mass transfer is due 


to a focused wind (ISidoli et al.ll2012ll . A focused 
wind or accretion stream requires a mass donor 
that nearly fills the Roche-lobe. If this is the case, 
we would expect variations in the mass accretion 
rate that would be attributable to a focuse d wind 
or accretion stream ( Blondin fc Owenlll997ll . This 
mechanism would lead to large variability in the 
X-ray intensity and could be observed in folded- 
light curves as well as hardness-intensity diagrams. 
Large intensity swings on the order of ~100 were 
observed in both Swift and X MM-Newton obser¬ 
vations of ^R J16418-4532 (|Sidoli et al.l 12012 : 
Drave et al. 2013ll . 


Near-infrared spectral features previously led 
to a spectral cl assification of eith er a late O- 
typ e supergiant ( Ghatv et al. 2008ll or a BN0.5 


la (Goleiro et al. 


20131 1. We note the presumed 
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Table 11 

Possible Parameters of Candidate Donor Stars for IGR J16418-4532 


Spectral Type 

M/Mq 

q 

R/Rq 

R'l/Rq^ 

Rh/R 

My 


E{J -Kf 

■Ob 

^sun 

(kpc) 

J d 

'-^sun 

(kpc) 

06 la 

44.10 

0.032 

19.95 

23.85 

0.84 

-6.38 

-0.21 

2.60 

77-81 

11.9 

11.9 

06.5 la 

41.20 

0.034 

20.22 

23.18 

0.87 

-6.38 

-0.21 

2.60 

73-76 

12.0 

12.1 

07 la 

38.44 

0.036 

20.49 

22.52 

0.91 

-6.38 

-0.21 

2.60 

70-73 

12.0 

12.3 

07.5 la 

36.00 

0.039 

20.79 

21.90 

0.95 

-6.38 

-0.21 

2.60 

67-69 

12.0 

12.4 

08 la 

33.72 

0.042 

21.10 

21.31 

0.99 

-6.38 

-0.21 

2.60 

64-66 

11.9 

12.6 


Note. —Possible Parameters of Candidate Donor Stars. 

“ The value for (J — K)q was calculated using (J — H)q and {H — K)o published in Martins fc Pled ( 20061) . E{J — K) is 
found by subtracting {J — K)o from the observed J — K 

^ The range of inclination angles of the system consistent with the measured eclipse half-angle. 

° The distance the object is from the Sun using the distance modulus. 

The distance the object is from the Sun using the radius to distance ratio derived from spectral energy distribution. The 
radius to distance ratio is found to be 3.77x10“^^ ( Chatv et al. 20081) . 


sis 




Fig. 19.— The black curves show the predicted 
eclipse half angle as a function of inclination angle 
for stars with the candidate spectral types for IGR 
J16418-4532. The red and black dashed lines indi¬ 
cate the eclipse half angle and estimated error as 
measured by Swift BAT. We assume a neutron star 
mass of 1.4 Mq (top) and of mass 1.9 Mq (bottom) 
and typical masses and radii for the assumed com¬ 
panion spectral type (see Table [T^. The blue ver¬ 
tical dashed lines indicate the lower limit of the 
inclination angle. Inclinations to the left of these 
correspond to stars that overfill the Roche-lobe. 


radius of both spectral types overfills the Roche- 
lobe, which shows that the proposed spectral clas¬ 
sifications are incorrect. The spectral type of the 
mass donor places an additional constraint on the 
source distance. Under the assumption that the 
K-band magnitude (juk) and extinction in the V- 
band (Ay) are magnit udes 11.48 and 14.5 for an 


08.5 I classification, ( Rahoui et al. 20081) found 


the distance of the source to be ^13kpc. Convert¬ 
i ng Ay to Ak using Table 3 in iRieke fc Lebofskv 
(119851 ). we confirm the calc ulations for the dis- 


(l200Sh 


tance of IGR J16418-4532 in lR.ahoui et al. 
using the distance modulus. The distance of IGR 
J16418-4532 assuming the aforementioned spec¬ 
tral types are reported in Ta ble ITT] u si ng th e val- 
ues for My ob t ained from IWegnen (120061) and 
Martins fc Pled ( 20061) . We also use the radius 


to distance ratio from the near-infrared spectral 
ene rgy distribution (SE P) measurements reported 
in ( Rahoui et al. 20081) together with our eclipse 
measurements to estimate the distance of IGR 
J16418-4532. We find the distance to be between 
11.9-12.6k pc, which is sl i ghtly smaller than that 
reported in Rahoui et al. ( 20081) for an 08 I star. 

IGR J16418-4532 is a heavily absorbed SFXT 
where the observed measured to be between 
3.9j)n'o> < 10^^ atoms a nd 7±2x 10^^ atoms 

(jRomano et al.l 120121 ). exceeds the val- 


cm 
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ues r eported by the Leide n/Argentine/Bonn sur- 
vey ( Kalber]a.eL_a]J 2005) and in the review by 
Dickey fc Lockmanl ( 1990 . which are 1.59x10^^ 
and 1.88x10^^ atoms cm“^, respectively. Since 
the measured Nn was found to be in excess of the 
Galactic H I, determining the inters tellar fraction 
of iVn is problematic and the value in iR.ahoui et al.l 
(|2008|) for the extinction cannot be verified with¬ 
out the presence of systematic error. 

We compare the observed value of the J — K 
color of 2.3ljl with the intrinsic {J — K)q for the 
proposed mass donors (see Table HH). Calculat¬ 
ing the difference between the observed J — K 
and the intrinsic (J — we find the redden¬ 

ing values E[J — K) for each proposed spectral 
type for the mass donor (see Table fTTI) . We calcu¬ 
late the reddening i n the B — V band, EjB — V), 
using Equation 1 in iGiiver fc Oumlzell (l2009ll and 
the e xtinction in the V-band (Ay) ( Rahoui et al.l 
2008 1. Conver ting E(B — V) to EjJ — K) us¬ 


ing Table 3 in Rieke fc Lebofskvl (198^, we find 


E{J — K) to be 2.45. While this is slightly lower 
than what would be expected for late O supergiant 
spectral types (see Table [TT|l. the presence of sys¬ 
tematic error described above prevents an exact 
calculation. 

These results show that stars with spectral type 
07.5 I and earlier satisfy the constraints imposed 
by both the duration of the eclipse and the Roche- 
lobe (see Figure [H]). Since the measured A^h is 
lar gely in excess o f inte rstell ar values measured 


lar gely m excess o t inte rstell ar values measured 
by iKalberla et all ( 2005 1 and Dickev fc Lockman 
( 1990l l. determining what fraction of TVh is inter¬ 
stellar in origin is problematic. We find measure¬ 
ments such as the distance to be consistent with 


the distances determined bv iRahoui et al.l (|2008fl . 


Given our measurements, spectral types near 07.5 
I are reasonable. By constraining this type of star, 
we solved the first part of a three part problem. 
The remaining pieces include pulse-timing mea¬ 
surements which can be done in the X-ray and 
radial velocity measurements in the near-infrared. 

4.1.3. IGR J16479-4514 

Our results show that while the expected 
masses and radii for the previously proposed spec¬ 
tral types for IGR J16479-4514 satisfy the eclipse 
half-angle (see Table | 8 ]), the implied radius would 


be la rger than the Roche-lobe radius (|Eggleton 
198311 . This is similar to the situation observed 
in IGR J16418-4532, which we describe in Sec¬ 
tion HXH We calculate the Roche-lobe radius for 
stars of spectral types 08.5 I and 09.5 lab to be 
19.1 i ?0 and 18.11?©, respectively (see Table [ 8 ]). 
Since the radius for each proposed spectral type 
is too large to satisfy the constraint imposed by 
the Roche-lobe (see Table | 8 ]), we suggest that the 
donor must be a slightly earlier spectral type. 

We find that spectral classes 07 I and earlier 
satisfy the Roche-lobe constraint (see Table IT^ . 
The ratio between the radius of the donor star 
and that of the Roche-lobe, /3, is found to ex¬ 
ceed 0.9, which is cons istent with other HMXB s 
that host supergiants ( Joss fc Rappapor'^ 1984 1. 
Transitional Roche-lobe has been proposed in IGR 
J16479-4514 ( Sidoli et al.l 2013ll . which requires 


a mass donor th at nearly fills the Roche-lobe. 


Sidoli et al.l (|2013ll found phase-locked flares were 


found in their observation of IGR J16479-4514. 


Sidoli et al.l (|2013ll attributes the physical mecha¬ 


nism responsible for these flares which are spaced 
0.2 in phase as large-scale structures in the wind. 


Usi ng near-infrared spectral features. iGhatv et al 


( 2008ll previously determined the spectral classifi¬ 
cation ofthe_donorstOT to be a late 0 -type super¬ 
giant ( Chatv et al.ll2008ll . The donor spectral type 
places an additional constraint on the distance of 
IGR J16479-4514. The K-band magnitude (tok) 
and the extinction in the V-band (Ay) a re found 
to be 9.79 and 18.5 ( Rahoui et al.ll 200 ^ . Addi- 
tionally, the ??« / ratio was found to be 1 x 10 “^° 


(IGhatv et al.l 120081 1. Using this information, the 


minimum distan ce of IGR J16479-4 514 was found 
to be '^d.Okpc ( Ghatv et al. 2008ll. Convert ini 


Ay to Ak using Table 3 in iRieke fc Lebofsk 


(|l985ll . we find Ak to be 2.07 and confirm the 
calculation for the distance of IGR J16479-4514 
using the distance modulus. The distances of 
IGR J16479-4514 assuming the aforementioned 
spectral types are reported in Table m usin g the 
values for My obtaine d from Wegneil ( 2006h and 
Martins fc PlezI (200^. We find that the distance 


of stars with spectral type 07 I and earlier to 
be between 4.4-4.6kpc, which are r esonably con- 


sistent with the measurements in Rahoui et al 


^http://www.iasfbo.inaf.it/-masetti/IGR/sources/16418.html 


(1200811 . 

IGR J16479-4514 is a heavily absorbed SFXT 
where the AA was measured to be 9.5±0.3x 10^^ atoms 
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Table 12 

Possible Parameters of Candidate Donor Stars for IGR J16479-4514 


Spectral Type 

M/Mq 

q 

R/Rq 

R'l/Rq^ 

Rh/R 

Mv 


E{J -Kf 

■Ob 

^sun 

(kpc) 

J d 
'^sun 

(kpc) 

06 la 

44.10 

0.032 

19.95 

22.03 

0.91 

-6.38 

-0.21 

3.48 

69-75 

4.44 

4.50 

06.5 la 

41.20 

0.034 

20.22 

21.42 

0.94 

-6.38 

-0.21 

3.48 

66-71 

4.46 

4.56 

07 la 

38.44 

0.036 

20.49 

20.80 

0.99 

-6.38 

-0.21 

3.48 

63-68 

4.46 

4.62 

07.5 la 

36.00 

0.039 

20.79 

20.23 

1.03 

-6.38 

-0.21 

3.48 

60-64 

4.48 

4.69 

08 la 

33.72 

0.042 

21.10 

19.68 

1.07 

-6.36 

-0.21 

3.48 

58-61 

4.44 

4.81 


Note. —Possible Parameters of Candidate Donor Stars. 

“ The value for (J — K)q was calculated using (J — H)q and {H — K)o published in Martins fc Pled ( 20061) . E{J — K) is 
found by subtracting {J — K)o from the observed J — K 

^ The range of inclination angles of the system consistent with the measured eclipse half-angle. 

° The distance the object is from the Sun using the distance modulus. 

The distance the object is from the Sun using the radius to dista nce ratio derived from spectral energy distribution. The 
radius to distance ratio is found to be 1x10“^° ( Chatv et al. 20081) . 


cm ^ (|Sidoli et al. 


201,11 ). This is an order of 
magnitude larger than the Galactic H I val¬ 
ues r eported by the Leide n/Argentine/Bonn sur- 
vev (iKalberla et al.l 2005) and in the review by 
Dickey fc LockmanI ( 199C ). which are 1.87x10^^ 
and 2.14x10^^atoms cm“^, respectively. Since 
the measured iVn for both sources was found to 
be in excess of the Galactic H I, determining the 
interstellar fraction of iVn is difficult. Therefore, 
we cannot verify th e value for the extinction in 


Rahoui et al.l (120081 ) without the presence of sys¬ 


tematic error. 

We calculate the reddening values E{J — K) for 
each proposed spectral type for the mass donor 
(see Table [T2|) using the observed value of the J — 
K color of 3. 270 and the intrinsic {J — K)q for the 
proposed mass donors (see Table [I2]). To check 
for consistency with late O supergiant stars, we 
compare this with t he reddening in the B—V band 
using Equation 1 in Giiver fc Oumlzel (l2009l) and 
the e xtinction in the V-band l7l\^) (lR,ahoui et al 


20081 ). Goi iverting E(B — V) to E( J — K) using 
Table 3 in iRieke fc Lebofs^ ( 198,51 ). we find the 


reddening in the E{J — K) band to be 3.12. While 
this is slightly lower than what would be expected 
for late O supergiant spectral types, the presence 


of systematic errors, described above prevents an 
exact calculation. 

Stars with spectral type 07 I and earlier satisfy 
the eclipse duration and Roche-lobe constraints 
(see Figure 1^01) . Determining the interstellar frac¬ 
tion of TVh is difficult since the measured A^h is 
greatly in excess of interstellar values. We find 
the distance of the proposed count erparts to be 


consis tent with those determined bv iRahoui et al 


■^http: //www. iasf bo. inaf. it/^masetti/IGR/sources/16479. html 


(|2008D . Since no pulsation period was identified, 
the next step to constrain the donor star would be 
radial velocity measurements in the near-infrared. 

4.14. IGR J18027-2016 

Our results show that the mass and radius 
of the donor star in IGR J18027-2016 can be 
constrained to be between 18.6±O.9M0 and 
17.4±0.9 i?Q and 19.4±O.9M0 and 19.5t;(;?Ro 
for edge-on orbits and inclinations where the 
Roche-lobe is completely filled, respectively. We 
find the inclination where the Roche-lobe is filled 
to be 73.3°, wh ich is consis t ent w ith the ear¬ 
lier results from iMason et al.l (120111) . Since the 
semi-amplitude of the radial velocities of both 
the donor star and compact object are known, 
we also find that the mass of the neutron star 
can be constrained. We calculate the mass of the 
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Fig. 20.— The black curves show the predicted 
eclipse half angle as a function of inclination angle 
for stars with the candidate spectral types for IGR 
J16479-4514. The red and black dashed lines indi¬ 
cate the eclipse half angle and estimated error as 
measured by Swift BAT. We assume a neutron star 
mass of 1.4 Mq (top) and of mass 1.9 Mq (bottom) 
and typical masses and radii for the assumed com¬ 
panion spectral type (see Table [T^. The blue ver¬ 
tical dashed lines indicate the lower limit of the 
inclination angle. Inclinations to the left of these 
correspond to stars that overfill the Roche-lobe. 


neutron star to be between 1.37±O.19M0 and 
1.43±O.2OM0 for our lower and upper limits (see 
Section 13.41) . 


Since the radius of the donor star is constrained, 
we also can estimate the distance, optical and 
near infrared magnitudes of IGR J1 8027-2016. 


Using SED measurements, iGhatv et al.l ( 2008h 
calculated the radius to distance ratio to be 
4±lxl0“^^, where the uncertainties are at the 
90% confidence interval. At the 90% confidence 
interval, we find the eclipse half-angle to be 34^3° 
and the radius of the donor star to be constrained 
to between 17t.i 19±1 Rq in the two lim¬ 

its described in Section l3i4l Combining our results 
for the radius of the do nor star with the radius to 
distance ratio (Table 6: IChatv et al.llioosh . we find 
that the distance of IGR J18027-2016 can be con¬ 
strained to ll±2kpc and 12±2kpc i n the two lim¬ 


its. U sing the distance modulus fe.g. lNespoli et al 


2008ll . the absolute magnitude of IGR J18027-2016 
can be calculated. The apparent magnitude in 
the R-band, R, and the extinction in the V-band, 


Ay w ere found to be 16.9 and 8.5 (jMasetti et al 


20081 1. We find the absolute magnitude in the R- 
band, Mr to be ^-5 at both constraints, which is 
what would be expected for a B-typ e supergiant 
( Martins fc Pled 2006 : Wegner 2006 1. 

Our results a re consistent with th e previously 
proposed B1 Ib (Torreion et al. 20inll or BO-Bl I 
( Mason et al. 2011 1 spectral types in IGR J18027- 
2016. We constrain the mass and radius of the 
donor star to be between 18.6-19.4M 0 and 17.4- 
19.5 7?0. We also constrain the mass of the neu¬ 
tron star to be 1.18-1.63 M^n, which ma r ginall y 


improves on the results in iMason et al.l (|20111 1. 


which was found to be 1.15-1.85M 0 . 

4.1.5. XTE J1855-026 

We find that the mass and radius of the donor 
star in XTE J1855-026 are constrained between 
19.6±1.1M0 and 21.5±O.5i?0 at edge-on orbits to 
2O.2±1.2M0 and 23.O±O.5i?0 where the Roche- 
lobe is just filled. The inclination where the 
Roche-lobe is filled is found to be 76.4°. We 
find the derived mass es and radii to be consi stent 
with those re ported in Carroll fc Ostliel ( 2006h and 
Allen ( 2 OOOII for stars with spectral type BO I. 
Since the semi-amplitude of the radial velocities 
for both the donor star and compact object are 
known, we find that the mass of the neutron star 
can be constrained to be between 1.77±O.55M0 
and 1.82±0.57 Mq (see Section lXSl) . We note that 
the large error estimates in the mass of the neu¬ 
tron star are likely to be attributed to substantial 
un certainties in th e estim ate in Ko as reported 
Gonzalez-GalanI ( 2015 ). This is likely to be 


attributed to emission line contamination and/or 
changes in the stellar wind. 

Based on optical and near-infra red spectra to¬ 
gether with the analysis reported in IVerrecchia et al 
(I 2 OO 2 I ). the spectral type of the donor star 
was previously determ ined to be a BO laep 
( Negueruela et al.ll2008l) . Based on the ratio of 
the equivalent widths of Si IV to Si III which is a 
diagn ostic for supergiant spectral types (jWalborn 
19711 ). the mas s donor was rec e ntly r efined to 

a BN0.2 la by iGonzale z-Galan ( 2015 ). _ Using 

SED measurements, Coleiro fc Chatvl ( 20131) cal¬ 
culated the radius and distance to be 26.9 Rq and 
10.8±1.0kpc, which places XTE J1855-026 in the 
Scutum arm region. Combining the properties 
of the newly derived spectral type with the dis- 
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tance modulus, iGonzalez-Galad (120151) recently 
calculated the nominal distance of XTE J1855- 
026 to be 101*14 kpc. Based on these results, the 
radius to distance ratio can be calculated to be 
5.6±0.5xl0“^^. Combining our results for the ra¬ 
dius of the donor star with the calculated radius 
to distance ratio, we find the distance of XTE 
J1855-026 can be constrained to 8.6±0.8kpc and 
9.2±0.9kpc in the two limits. 

Corbet et al. ( 19991) measured the IVh in XTE 


J1855-026 to be 14.7±0.6x 10^^ atoms cm”*^, 
which exceeds the values r eported by the Lei- 
den/Argentine/Bonn survey (Kalberla et al.l2005 ) 
and in the review by IPickev fc LockmanI ( 1990f ). 
In a Swift XRT observation, the A^h was found 
to be 4.1±0.5x10^^atoms cm“^ (Romano et al 


20081 ) ■ These are signifcantly larger than the 
measurements for the i nterstellar AAr, which are 


6.62x10^^ atoms cm ^ (Kalberla et al. 2005) and 


7.35x10^^ atoms cm ^ ( Dickev fc Lockmanlll99^ . 


Since the measured A^h was found to be in excess 
of the Galactic H I, the conversion between the 
interstellar A^h and the extinctio n Ay is problem¬ 
atic and the value of 5.8±0.9 in IColeiro fc Ghatv 


()2013l) cannot be verified without the presence of 
systematic error. 

4.2. What is the nature of the non-zero 
eclipse flux in IGR J16393-4643? 

The source flux during eclipse in IGR J16393- 
4643 does not reach 0 counts s“^ in the folded light 
curves (see Figure |4]). We find the ratio between 
the flux during eclipse to the flux outside eclipse to 
be 54±5%. This is significantly larger than what 
is observed in the other XRBs in our study (see 
Table HI). 

We first discuss the possible scenario where we 
observe a partial eclipse in IGR J16393-4643. In 
our model for the eclipse (see Section IQ) , we 
assume the compact obj ect to be a point source 
(jjoss fc Rappaportl Il984l) . While this is a valid 
assumption since the radius of the compact object 
is much smaller than that of the donor star, our 
approximation does not consider the extended X- 
ray emission region. In this case, the constraints 
the mass and radius of the donor star must take 
into account the size of the extended emission re¬ 
gion. Since a significant residual flux is observed, 
it is likely that we observe a partial eclipse. 


We also consider the possibility that the resid¬ 
ual emission observed in IGR J16393-4643 is 
attributed to a dust-scattering halo similar to 
what is observed in some othe r HMXBs (Gen 
X-3; Vela X-1; GAP 1657-415. IPav fc Tennant 


19911 : IWoo et al.lll994lAudlev et al.ll2006l) . While 


residual emisison has been attributed to a dust¬ 
scattering in some other HMXBs, the residual 
emission found in IGR J16393-4643 in the BAT 
folded light curves (15-50 keV) is seen at much 
higher levels (see Tabled]). A dust-scattering halo 
is predominantly a soft X-ray phenomenon. While 
a significant fraction of the out -of-eclipse flux may 
be in a dust-scattering halo (|Predehl fc Schmitt 
Il99,5l) . we expect a smaller fra ction to be in the 
energies resolved by BAT (e.g. iGolev et al.l 12014 , 
and references therein). Therefore, we conclude a 
dust-scattering halo cannot be the sole mechanism 
responsible for the residual emission seen in the 
folded light curves. 

Finally, we discuss the possibility that Comp¬ 
ton scattering and reprocessing in a region of dense 
gas could accou n t for t he residual flux in eclipse. 


Bodaghee et al.l ()2006l ) found that a Compton 


emission (comptt within Xspec) model with an 
electron temperature of 4.4±0.3keV and optical 
depth of 9±1 provides a good fit to the average 
spectrum of IGR J16393-4643. Additionally, Fe 
Ka and Fe K/3 lines were found at 6.4 keV and 
7.1keV, respectively, where the ratio of the iron 
intensiti es is seen to be consist e nt with photoion¬ 


ization ( Bodaghee et al. 2006t Kaastra fc Mewe 


1 19931 : llslam et al. 2015h . The equivalent widths 


(EQW) in the XMM-Newton observation were 
found to be 60±30eV and an upper limit of 
120eV . for the Fe Kg and Fe K/3 lines respec¬ 
tively (|Bodaghee et al.ll2006l) . In a recent Suzaku 
observation, the EQW of the Fe Ka and Fe K/3 
lines wer e found to be 46^ 1 eV and an upper limit 
of30eV (jlslam et al.ll2ni5[) . While the mechanism 
responsible for the Fe Ka and Fe K/? is likely to be 
fluorescence of cold matter, the equivalent widths 
point to a like ly origin in a spherica l distribution 


of dense gas ( Bodaghee et al. 20061) . Therefore 


Compton scattering and reprocessing might be 
the sole contributor to the residual emission ob¬ 
served in the BAT folded light curves. 

It is likely that a combination of the mech¬ 
anisms described above account for the residual 
emission found in the folded light curves of IGR 
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J16393-4643 (see Section 13.11) . Since the count 
rate during eclipse is significantly larger than what 
is observed in most other eclipsing HMXBs, it is 
likely that only a small fraction comes from a dust¬ 
scattering halo. 


4.3. What mechanism is responsible for 
asymmetries in the eclipse profile? 


The eclipse profiles show the presence of asym- 
metries (see T a bles[^ and[3t as previously noted by 


Falanga et al.l ()2015|) in the cas es of IGR J18027 - 


2016 and XTE J1855-026 and I.Tain et al.l (|2009^ 


in the case of IGR J16479-4514. These asymme¬ 
tries seen in the ingress and egress durations are 
suggestive of the presence of complex structures 
in the wind such as accretion or photoionization 
wakes. 


We first discuss the possible case that the asym¬ 
metry in the eclipse profiles can be attributed to 
accretion wakes. In an HMXB driven entirely by 
a spherical wind, material is only accreted in a 
cylindrical region where the kinetic energy is less 
than the gravitational potential energy of the com¬ 
pact object. The radius of the accretion cylin¬ 
der is the Bondi-Hoyle accr etion radius (Equa¬ 
tion 1, iFeldmeier et al. 1996ll . Perturbed mate¬ 
rial forms an “accretion wake” that typically trails 
the orbit of the neutron sta r and results in larg e 
intrinsic column densities ( Blondin et al 
Prior to eclipse, the progressively increasing A^h 
partially obscures the X-ray emission resulting in 
longer ingress durations compared to the dura¬ 
tion of egress. Since the accretion wake is located 
beyond the compact object during egress, no ap¬ 
parent increase in the intrinsic A^h is observed. 
The ingresses observed are somewhat larger than 
egress, which is consistent with the p resence of an 


accretion wake (Blondin et al. 199011. The count 


rate prior to ingress is also somewhat smaller than 
that observed after egress, providing additional 
evidence for accretion wakes. Hardness ratios or 
measurements of A^h folded on the orbital period 
could be implemented to confirm the possibility of 
accretion wakes. 


We also consider the possibility that photoion¬ 
ization wakes could explain the asymmetric eclipse 
profiles. The eclipse profiles of IGR J18027-2016 
and XTE J1855-026 are compared to those seen 
in eclipsing systems where asymmetric density en¬ 
hancements are observed on large spatial dimen¬ 


sions (e.g. Vela X-1: IFeldmeier et al.lll996ll . In 
systems where the X-ray luminosity is significantly 
high, a switch-off of the radiative driving force 
could lead to a red uced wind velocity a nd en¬ 
hanced wind density ( Feldmeier et al.iri996[) . This 
enhanced X-ray scattering region trails the neu¬ 
tron star and results in ingress durations that 
are sig nificantly larger than those observed at 


egress (jFeldmeier et al.Ml996[ l. The eclipse pro¬ 


files of IGR J18027-2016 and XTE J1855-026 differ 
from those expected from a photoionization wake- 
the ingre ss duration in Vela X -1 was seen to be 
(j) =0.11 ( Feldmeier et al.l 19961) . 

We additionally discuss how energy dependence 
in the asymmetric eclipse profiles can arise. The 
high Vh, on the order of 10^^ atoms cm“^, im- 
plies the presence of a strong circumstellar wind 


Hr K 

(lKuulkersll2005ll , and the X-ray absorption due to 
this will lead to sharper ingress and egress tran¬ 
sitions compared to that seen at lower energies. 
In a Suzaku observation of IGR J16479-4514 that 
covers part of the eclipse, Sidoli et al. ( 2013ll found 
the egress trans ition to be broader in soft energies. 
In their review. IPalanga et al. ( 2015[) found asym¬ 
metries to be slightly enhanced at lower energies 
compared to higher energies. While investigating 
the energy dependence in the eclipse transitions is 
beyond the scope of the present paper, this study 
in obscured SGXBs will be difficult due to a re¬ 
duced count rate at low energies that results from 
the la rge intrinsic absorpt ion present in these sys¬ 
tems ( Falanga et al. 20151) . 

Finally, we consider the possibility that the 
asymmetric eclipse profiles could be attributed 
to a small to modest eccentricity. Since we con¬ 
sidered objects with relatively short orbital peri¬ 
ods, we expe ct the eccentricity of the systems t o 
be near zero ( Zahnl[l97^ iMaccarone et al. 2014h . 
The eccentricity in IGR J18027-2016 and XTE 
JI855-026 were both noted to be small to modest, 
where e was found to be less than 0.2 and 0.04, 
respe ctively (|Augello et al.ll2003t iGorbet fc Mukai 
2 OO 2 II . In the cases of IGR J16393-4643, IGR 
JI6418-4532 and IGR JI6479-4514 where no 
pulse-timing or radial-velocity methods are avail¬ 
able to determine an orbital solution, we constrain 
the maximum allowed eccentricity to that where 
the radius of the donor st ar completely fi l ls the 


Roche-lobe at periastron (jGoossens et al.l 12013 . 


and references therein). These were all found to 
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be near zero (see Figure [2T]) . Since the eccentric¬ 
ities are found to be small to modest, we do not 
expect these to result in sizeable asymmetries in 
the eclipse profile. Additionally, apsidal advance 
will be apparent in the case that an eccentric orbit 
could lead to asymmetries in the eclipse profile. 
While apsidal advance would be difficult to detect 
in the ^9yr of Swift data, we believe it to be un¬ 
likely for asymmetries to be solely attributed to 
small to modest eccentricities. Furthermore, accu¬ 
rate measurements of apsidal advance will depend 
on multiple pulse-timing measurements of these 
systems, which are not yet available. 



Phase 


Fig. 21.— A sample plot of the Roche-lobe of 
the donor stars in IGR J16393-4643 (top), IGR 
J16418-4532 (middle) and IGR J16479-4514 (bot¬ 
tom) vs. orbital phase for eccentricities ranging 
from 0.0 to 0.1 in steps of 0.02. The horizontal red 
line represents the radius of the donor star under 
the assumption that the donor is an 09 I for IGR 


J16393-4643 and 07 I for both IGR J1641 8-4532 
and IGR J16479-4532 (|Martins et al.ll2005ll . 


5. Conclusion 


Eclipsing X-ray binaries provide an opportunity 
to constrain the physical parameters of the donor 
star as well as the compact object. To determine 
the eclipse half-angle in our survey, we modeled 
the eclipses using both symmetric and asymmet¬ 
ric step-and-ramp functions. The luminosity of 
each syst em is less than ex pected for Roche-lobe 
overflow ( Kaner et al. 2004^ , which means we can 
attach the constraint that the mass donor under¬ 
fills the Roche-lobe. Since IGR J18027-2016 and 


XTE J1855-026 are the only “double-lined bina¬ 
ries” in our sample, we calculate the parameters of 
the other systems assuming the neutron stars to be 
at the white-dwarf Ghandrasekhar limit, 1.4M0. 
We also calculated the parameters of the other 
systems assuming a more massive neutron star- 


1.9M©. 

Our results show that stars with spectral type B 
III satisfy both constraints imposed by the eclipse 
duration and Roche-lobe for IGR J16393-4643. 
Assuming the estimates f or the mass and r a dius o f 
a B5 III star reported in Garroll fc Ostliel (2006), 
we find the mass and radius of the donor star to 
exceed 7 Mq and 6.3 Rq. B I stars were found 
to overfill the Roche-lobe. The source emission 
in IGR J16393-4643 does not reach 0 counts s“^ 
in the folded light curves, where the fraction be¬ 
tween the flux in eclipse to that outside eclipse was 
found to be 54±5% (see Tables [T]-13|) . Compton 
scattering and reprocessing in a dense region of 
gas could possibly account for the X-ray emission 
region not obscured by the donor star. 


Our results show that the previously proposed 
08.5 I and BN0.5 la spectral types for the mass 
donor in IGR J16418-4532 must be excluded. 
While these spectral types satisfy the eclipse half¬ 
angle, the Roche-lobe is significantly overfilled. 
Stars with spectral type 07.5 I or earlier are con¬ 
sistent with both the eclipse half-angle and the 
Roche-lobe. In this case, we find the mass and 
radius of the donor star to exceed 36.00 M© and 
20.79 i?© assuming the estimates fo r the mass and 
radius of an 07.5 I star reported in iMartins et al.l 
(120051 1 . We find the minimum inclination angle 
of the system to be 67°. The distance measure¬ 
ments of IGR JI64I8-4532 are consistent with the 
previously determined distance (see Table HD); 
however, determining the interstellar fraction of 
A^h was found to be problematic. The ingress 
and egress durations in the folded light curves 
were found to be asymmetric where the ingress 
duration is longer than the egress duration. This 
is likely attributed to the presence of an accre- 
tion wake or from a f ocused stream as n oted in 
Sidoli et "ni ( 20I2ll and Drave et al. ( 2013 1. 


The previously proposed 08.5 la and 09.5 lab 
spectral classifications of the mass donor in IGR 
J16479-4514 must be excluded because the Roche- 
lobe is significantly overfilled. However, we found 
that stars with spectral type 071 and earlier sat- 
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isfy both constraints imposed by the eclipse du¬ 
ration and Roche-lobe. Assuming the estimates 
fo r the mass and radius of an 07 I star reported 
(l2005ll . the mass and radius of 


Martins et al 


the donor star are found to exceed 38.44 Mq and 
20.49i? 0 , respectively. We find the minimum in¬ 
clination angle of the system to be 63°. The dis¬ 
tance measurements remain unchanged from ear¬ 
lier measurements (see Table [T2|); however, the in¬ 
terstellar fraction of A^h remains undetermined. 
We find the ingress and egress durations to be 
symmetric within the error bars. The ratio be¬ 
tween the radius of the donor star and Roche-lobe 
was found to exceed 0.9, which shows the possibil¬ 
ity of transitional Roche-lobe overflow. 


The mass and radius of the donor star in IGR 
J18027-2016 was constrained to 18.6±O.9M0 and 
17.4±O.9i?0 and 19.4±O.9M0 and 19.5t!^;?R0 
in the two limits. We find the inclination angle 
where the donor star just fills the Roche-lobe size 
to be 73.3°. We also find the distance measure¬ 
ments of IGR J18027-2016 to be ll±2kpc and 
12±2kpc in the allowed limits. In the allowed 
limits, we constrained mass of the neutron star to 
between 1 . 37 ±O.I 9 M 0 and 1 . 43 ±O. 2 OM 0 . The 
folded light curve shows complicated and asym¬ 
metric ingress and egress durations, which can 
be explained either by the presence of accretion 
wakes. 


Our results show the mass and radius of the 
donor star in XTE J1855-026 to be constrained to 
19.6±I.I Mq and 2I.5±0.5 Rq at edge-on orbits to 
2 O. 2 ±I. 2 M 0 and 23.O±O.5i?0 where the Roche- 
lobe size is just filled. We find the inclination an¬ 
gle where the donor star just fills the Roche-lobe 
size to be 76.4 degrees. In the allowed limits, we 
find the distance of XTE JI855-026 can be con¬ 
strained to 8.6±0.8kpc and 9.2±0.9kpc. We find 
the mass of the neutron to be constrained between 
1 . 77 TO. 55 M 0 and 1.82±O.57M0. Complicated 
and asymmetric ingress and egress durations were 
seen in the folded light curve, which suggests the 
presence of complex structure in the wind. 

To further constrain the physical parameters 
of the donor star and the compact object in all 
these systems, additional observations are re¬ 
quired. Constraining the mass of the neutron star 
will help constrain the Equation-of-State. Since 
the pulse period has been accurately measured, 
for IGR JI6393-4643 and IGR JI64I8-4532 the 


study would benefit from both pulse-timing anal¬ 
ysis as well as radial-velocity curves in the near- 
infrared. A radial-velocity curve in the optical or 
near-infrared would provide additional constraints 
for IGR JI6479-4514 where a pulse-period has yet 
to be identified. 

We thank the anonymous referee for useful com¬ 
ments and support from NASA 14-ADAP14-0I67. 
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